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ON  THE  USE  OF  FINITE  FIELD  BASED  MODELING 
IN  POLYNOMIAL  MANIPULATION 

Stephen  John  Nuspl,  Ph0D» 
Department  of  Electrical  Engineering 
University  of  Illinois,  1968 

A  major  hindrance  to  the  manipulation  of  polynomials 
by  a  computer  is  the  rapid  growth  of  the  polynomial  coeffi- 
cients during  the  application  of  some  frequently  used  algor- 
ithms o   The  methods  presented  in  this  thesis  to  try  to  alle- 
viate this  problem  and  to  increase  the  efficiency  of  polyno- 
mial manipulation  are  based  on  the  use  of  algebraic  redundancy 

Let  Iq  be  the  set  of  all  integers  and  I«+i   d© 
the  adjunction  of  I.   consisting  of  the  set  of  all  polyno- 
mials in  x*+1  with  coefficients  in  loo  Let  FQ  be  a 
field  of  order  pQ  where  pQ  is  a  prime  in  IQo  The  exten- 
sion field  P,+1  is  defined  to  be  the  set  of  residues  in 
F.[x.+,]  modulo  p.  ,„   a  monic  irreducible  polynomial  in 
^1^x1+l^°   ^he  sets  I*   are  defined  as  the  modeled  domains 
and  the  sets  F,  and  their  adjunctions  as  the  model  domains* 
A  modeling  function  which  maps  any  element  of  a  modeled 
domain  into  any  one  of  the  model  domains  in  defined  and 
shown  to  be  a  homomorphisnu 

Some  of  the  models  are  capable  of  retaining  the 
degree  information  of  selected  indeterminates  in  the  modeled 
polynomials o  A  virtual  degree  augmentation  is  added  to  the 


model  and  its  algebraic  treatment  is  developed  in  such  a  way 
that  it  is  able  to  detect  when  a  model  polynomial  has  the 
same  degree  as  the  modeled  polynomial  without  looking  at  the 
latter o   A  major  result  based  on  the  augmented  modeling  is 
the  model  conclusive  g.c.d.  theorem  which  gives  the  conditions 
under  which  two  polynomials  can  be  shown  relatively  prime  by 
performing  the  g0Cod»  algorithm  only  on  their  models »   Since 
model  domains  in  which  the  coefficients  are  finite  field 
elements  can  be  chosen,  the  computational  effort  required 
can  in  many  cases  be  orders  of  magnitude  less  than  if  the 
goC.do  algorithm  were  applied  directly  to  the  modeled 
polynomials o  Further  applications  of  modeling  are  discussed 
but  not  developed  to  any  depth* 
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CHAPTER  I   INTRODUCTION 

The  increasing  use  of  the  computer  as  an  aid  in 
algebraic  manipulation  makes  it  desirable  to  study  techniques 
by  which  the  speed  of  performing  a  manipulation  is  increased 
and  the  amount  of  physical  equipment  required  during  the  mani- 
pulation is  decreased o   Of  central  importance  are  techniques 
by  which  polynomials  in  several  variables  are  handled  since 
most  problems  eventually  filter  down  to  a  total  description 
by  polynomials  or  a  description  in  which  polynomials  play  a 
key  roleo   In  most  practical  cases  the  polynomials  are  defined 
over  the  integral  domain  of  integers  and  consequently  it  is 
algorithms  in  these  polynomial  domains  which  should  be 
optimized  for  speed  and  minimality  of  equipment o 

The  method  proposed  in  this  thesis  to  accomplish  this 
is  to  add  to  the  internal  representation  of  the  polynomial 
redundant  information  in  such  a  way  that  the  redundancy  may 
be  used  to  determine  some  desired  properties  of  a  polynomial 
without  requiring  that  the  entire  polynomial  be  scanned  for 
the  property o   For  example „  we  may  have  two  polynomials  in  a 
number  of  variables  and  we  may  wish  to  ask  the  question  "Is 
polynomial  8a9  identical  to  polynomial  °b"  in  the  sense  that 
if  both  are  put  into  the  same  canonical  form  and  ordering 0 
for  every  term  in  "a1  there  is  an  identical  term  in  8b»?98 
The  straight  forward  way  of  answering  this  question  is  to 
put  "a"  and  "b8  into  the  same  canonical  form  and  compare  them 


term~by~term  until  a  mismatch  occurs  or  until  the  terms  are 
exhaust edo   When  such  an  algorithm  has  a  very  high  probability 
of  producing  the  answer  "no"  and  if  the  much  more  concise 
redundancy  is  able  to  state  conclusively  that  "no11  is  the 
answer,  considerable  total  time  can  be  saved  by  first  trying 
to  answer  the  question  only  with  the  redundancy „   Then  the 
actual  polynomials  must  be  used  only  if  the  redundancy  does 
not  produce  a  conclusive  "no" „   In  either  case  the  answer 
will  be  the  same  as  when  only  the  actual  polynomials  are  used 
but  because  of  the  high  probability  of  a  "no"  the  comparative 
average  time  for  answering  the  question  can  be  orders  of 
magnitude  lower  when  the  redundancy  is  used. 

Of  particular  interest  is  redundant  information 
which  retains  some  of  the  structure  of  the  polynomials „   When 
operations  are  performed  on  the  redundancy  and  the  result 
gives  some  information  about  the  effect  of  applying  the  equi- 
valent operations  on  the  actual  polynomials,  the  redundancy 
can  be  thought  of  as  being  a  "model"  of  the  actual  polynomial . 
Since  this  will  be  the  case  in  the  development  of  this  thesis „ 
we  will  define  the  sets  of  all  polynomials  over  the  integral 
domain  of  all  integers  as  the  "modeled  domains",  the  sets  of 
all  models  as  the  "model  domains"  and  the  process  of  mapping 
a  modeled  polynomial  into  its  model  as  a  "modeling  function90 « 

The  modeling  used  in  this  thesis  is  based  on  finite 
fields  and  their  adjunctions  in  the  sense  that  the  model 


domains  are  either  finite  fields  or  adjunctions  of  finite 
fields o   The  usefulness  of  this  particular  type  of  modeling 
stems  primarily  from  the  fact  that  the  modeling  functions 
which  will  be  defined  are  homomorphisms  from  the  modeled 
domains  into  the  model  domains  which  produces  the  desired 
effect  of  preserving  some  of  the  structure  during  the 
modeling o   It  will  be  shown  that  some  of  the  modelings  which 
will  be  defined  have  the  additional  property  of  being  capable 
of  retaining  the  degree  information  of  selected  variables 
and  burying  the  rest  in  the  finite  field  coefficients  of 
these  variables,,   The  practical  significance  of  this  is  that 
the  degree  information  of  the  selected  variables  can  be  used 
for  decisions  in  algorithms  while  at  the  same  time  the 
coefficients  are  finite  field  elements  which  Implies  that 
they  are  automatically  restricted  in  size0   Normally  the 
problem  in  polynomial  manipulation  is  that  the  coefficients 
grow  so  rapidly  during  the  application  of  an  algorithm  that 
it  is  impractical  to  work  with  other  than  almost  trivial 
polynomials  in  a  few  variables  even  when  a  very  fast  computer 
is  used  to  perform  the  manipulation,, 

The  only  other  effort  known  to  the  author  in  which 
finite  field  theory  was  used  as  an  aid  in  algebraic  manipula- 
tion is  ref0  (8)  by  W„  Ao  Martin0   On  the  one  hand  his  approach 
was  much  more  ambitious  than  that  contained  In  this  thesis 
since  the  modeling  was  extended  to  functions  of  a  complex 


variable o   However „  his  models  were  restricted  to  take  on 
values  only  in  a  prime  field  which  resulted  in  the  loss  of  all 
degree  information  in  the  models .   It  is  precisely  the  pre- 
servation of  this  information  which  gives  so  much  power  to  the 
modeling  developed  here  even  though  it  is  restricted  to  poly- 
nomials over  integers o   Finite  field  theory  has  been  found 
useful  in  determining  the  existence  of  solutions  of  equations , 
so  it  is  not  too  surprising  that  the  same  theory  can  also  play 
an  important  role  in  practical  manipulative  work  with 
polynomials o 

The  development  of  the  modeling  is  first  started  in 
Chapter  2  with  the  review  of  some  concepts  in  modern  algebra 
which  are  later  required 9  a  rigorous  definition  of  the  modeled 
and  model  domains  and  some  relations  among  these  domains »   The 
modeling  functions  are  defined  and  proven  homomorphisms  in 
Chapter  3°   The  detailed  homomorphism  proofs  are  included 
mainly  for  completeness s  they  need  not  be  studied  in  order  to 
understand  the  modelings 0   Finally,  a  summary  which  should 
help  to  unify  and  to  make  more  understandable  the  relations 
between  the  various  modeling  functions  is  presented  in  section 
3o7o 

Chapter  4  introduces  the  additional  concept  of  a 
virtual  degree  augmentation  to  the  models „   A  set  of  operations 
defined  on  the  augmented  domains  are  developed  in  sufficient 
depth  to  be  able  to  arrive  at  an  adequate  set  of  functions 


which  permit  us  to  define  a  greatest  common  divisor  (g0Cod0) 
algorithm  which  operates  on  elements  in  the  model  domains <> 
It  is  then  shown  by  the  "model  conclusive  goCodo  theorem'8 
that  two  polynomials  can  be  shown  to  be  relatively  prime 
under  specified  conditions  by  using  only  elements  in  the 
model  domains  during  the  computation,.   The  practical  signi<= 
ficance  of  this  is  demonstrated  with  a  set  of  examples „ 

In  Chapter  5  some  practical  aspects  of  operations 
in  the  prime  subfield  and  in  other  model  domains  are  briefly 
considered o   The  appendix  includes  an  example  of  a  complete 
set  of  models  for  a  polynomial  in  3  variables  and  a  tabula= 
tion  of  a  set  of  primes  with  selected  properties „ 


CHAPTER  II   BACKGROUND,  POLYNOMIAL  DEFINITIONS,  DOMAINS 

2d  Introduction 

In  order  to  define  the  polynomial  modeling  process 
rigorously  it  is  first  necessary  to  review  the  definition  of 
a  polynomial  ring  and  the  operations  which  are  related  to  the 
modeling o   The  required  theory  Is  based  on  concepts  in  modern 
algebra  some  of  which  will  be  stated  in  section  2o2»   Much  of 
the  material  is  very  elementary  for  someone  well  acquainted 
with  algebra o  but  its  inclusion  is  necessary  in  some  cases  for 
the  explicit  referencing  of  well<°known  theorems  and  In  other 
cases  it  Is  used  as  a  means  of  defining  the  notation  intrin- 
sically o  For  example 0    the  definition  of  a  group  Includes  the 
normally  used  existential  quantifier  (a)„   the  universal 
quantifier   (V)   and  the  symbol   (030)   standing  for  "such 
that"o   The  reader  not  familiar  with  a  concept  or  symbol  used 
should  consult  one  of  the  references  (1),  (10) „  (11) „ 

2o2  Notation,,  Conventions „  Polynomial  Rings 

Let   {S0O}   be  a  monoid  in  which  S  is  a  set  of 
elements  closed  under  the  binary  associative  operation 
O:  SXS  — >  S*   Let  0  and  1  represent  the  Identity 
element  In  S  when  the  operation  O   is  respectively  addi- 
tive and  multiplicative o  The  sum  symbol ,  Yo  will  then  be 
defined  as  follows? 


0  n+1         n 

For  a±   €  S0  2^0  (ai)  ■  aQ  g  Jo  (a±)   -   (  ^O  (a^Oa^^ 
i=0  1=0        1=0 

for  n  >  1 


When  O  Is  additive  we  make  the  additional  convention?  for 


a  a0  =  0 

a  €  S9  n  >  0  define  na  «  }o  (a.  )  where 

1^0  a4  =  a  0  <  i 

When  O  is  multiplicative  we  define 


an  -     )o  (a.  )  where   < 

ifro  \a1  =  a  0  <  i  . 

It  is  immediately  apparent  that  for  m  >  09  n  >  0 

ma  O  na  =  (m+n)a    when  O  is  additive 

amO  an  =  am+n      when  O  is  multiplicative, 


{O0O}      Is  a  group  if  G  is  a  non=>empty  set 

<2o2<>01> 

and  O  is  a  binary  operation  on  the  cartesian 

product  GXG  with  the  following  properties? 

1.  O  i  GXG  — >  G  (Closure) 

2*  a0b0©  €  G  s>  aO(bOc)  ■  (aOb)Oo-(Associatlvlty) 

3»  3e€Go9olfa€GcaOeaa 

^o  va€G03b€GoB8aObse 


8 

Definitions   {R0©0©}   is  a  ring  if   {R,©}   is  a  commutative 
<2o2o03>     group  and  ©  is  a  binary  operation  on  RX  R 
with  the  properties? 

lo  ©  §  RXR  >  R 

2o   a0b0c  e   R  =>  a  ©  (b  0  c)  »  (a  ©  b)  ©  c 
3o   a0b0c  6  R  ■>  a  ©  (b  ®  e)  ■  (a  ©  b)  ®  (a  ©  c) 
and   (a  ®  b)  ©  c  ■  (a  ©  c)  ©  (b  ©  c) 

The  ring   {R,®,©)  will  be  referred  to  as  "the  ring  R  with 
respect  to  the  operations  @0   ©w   or0  when  it  is  understood 
that  R  forms  a  ring  with  respect  to  the  operations  ®0  © 
simply  as  "the  ring  R"   or  "R"o 


Definitions   Let   {R  „■+„<>}„  {R10®0©}   be  rings  and 
<2o2o0^>     x    „  R  __>     R    Then  x  ls  a  rlng  nomo„ 

morphism  from  R  into  R,   if  ^a0b  6  R0  aX  ©  bX  =  (a  -t  b)X 
aX  ©  bX  =  (a  °  b)XQ  Under  these  conditions  we  denote  X 
by  X  s  R  -§->  Rxo 


Definition?   e  is  called  the  multiplicative  identity  or 
unit  of  R  if  e©a^a©e  =  a  for  each 
a  in  Ro   If  e  is  in  R9   R  is  said  to  be  a  ring  with 
an  identity o 


zt  =  {i  I  1   Is  an  integer  and  0  <  1} 


Definitions  An  infinite  R-sequence  is  defined  by  the  fol- 
<2o2o0?5     lowings  Let   {R0@9©}  be  a  ring  with  identity, 
Let  09  1  be  the  zero  and  unit  in  R0  An 


R=sequenee  is  then  defined  to  be  a  function  f  s  ZQ  *=— >  R 

where   (i)f  ■  t±     with  I  €  zj  t  ^  6  H,   f  is  denoted  by 

f  ■  (tqDt^0t^9 oo e tTs  9 o o o ) o  r.   is  sometimes  called  the 
i'th  coefficient  or  1'th  component  of  fQ 


Definition?  R^  «  [u  |  u  is  an  infinite  Resequence? 
<2c2o08>     u  m   (U0o11iOGOO)  ^a  s  n  «3.  ?  i  >  n,  ut  s  0} 

The  last  condition  is  equivalent  to  saying  that  only  a 
finite  number  of  the  components  of  u     are  not  the  zero 
element  in  Rc 


Definitions  Let  uev  €  H™  us  (u~,,u,9u 

<2o2o09> 


<F*   0  O  O    O    II    S 

c 


tJ  ui»vi  €  H°  We  defin« 


.». 


the  binary  operators  ©  f  ©  as  follows 


g  s  R^XR   ==>  R  where  u|fs  CuQ  ©  bo*-ul  ®  vl° 

®  %   R^Xr   — »>  R  where  u  ©  ▼  =  ^eo»Gl»G29  «  0  o ) 

i 

and  eA  s  £@  (^  ©  vi~i^  s   7  ®  ^uj  ®  vfc^ 


Theorems   {R-,*©,©}   is  a  ring  with  identity  in  which  the  zero 
<2o2ol0>  element  is   o  =  (0,0,0,,,,)   and  the  unit  is 


1  =  (1,0,0 


9000/0 


Definitions   R  =  {u  |  u  =  (u0,090,  .  0  8 )  s  uQ  €  R0  u  £  R-jJ 
<2o2oll> 


Theorems  R  c  R±     and   {lo®*®}  =  {R8@,®}  where  =  denotes 

<2  2  12> 
«.o«.o  t.   ring  isomorphism  and  with  the  isomorphism  defined 

by  the  function  X  s  R  — — >  R  where 

a  €  R  ->  aX  -   (a,0,0D,,J,   The  isomorphism  X     will  some- 
times be  referred  to  as  an  embedding  of  R  in  R,<, 

Definitions   Scalar  multiplication  is  defined  to  be  the 
<2o2ol3> 


binary  operation  9  *  RXR1  — — >  R1  where 
s>  a  ©  u  =  (a  ©  uQOa  ©  u^ 


a  6  Ro  u  6  R1   ~>  a  ©  u  =  (a  ©  uft,a  ©  u,  0  0 » o ) » 


Theorems   The  complex   {{R0®9©}9  {R,,®},  ©}   is  a  vector 

<?  ?  li^>  0   1? 

space  over  R  and   {x  ,x «>x  , «, » » }   is  a  set  of 

basis  elements  when  x  ■  (09100000 . » 0 )  €  R-,0 

Proofs   The  vector  space  properties  are  readily  proven?   of 

0  1 
interest  to  us  is  the  fact  that   {x  ,x  0OOO}      is  a  set  of 


basis  elements o   By  our  definition  of  exponentiation, 
.0 

i=0  i=0 


=  £  ©  (x)  ■  Is  x*3  =  V  ©  (x)  =  (0,0,oo,,l,0,o,o) 
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where  the  1  is  In  the   j°th  component  position  of  x' 
For  u  €  R-j^  u  =  {u0(u10u2„  c  . » )  0  ^  e  R9  u  can  be  ex~ 
pressed  as  u  =  )   @  (u*  ©  x  )•   It  follows  that 

0   1 

{x  „x  9ooo}   Is  a  set  of  basis  elements  for  R1 <> 


Definitions   If  s  €  R-.   and  a  ©  s  =  s  ©  a   fa  €  R.  s 

°  °  ->      is  said  to  be  a  scalar  over  Ro   Since  R 

is  commutative  every  element  in  R,   is  a  scalar  and 

of  particular  importance,,   x  is  a  scalar „ 


Definitions   v  ^s   said  to  be  algebraic  over  R  if 

v  €  R-j   and   In,  a,  0  0  <  i  <  n  such  that 

at  least  some  a^^  +  °  and  that   T  @  (a±  ©  v1}  =  0  , 

1^0 


Definitions  v  Is  said  to  be  an  indeterminate  over  R 
if  v  6  R-j,   and  v  is  not  algebraic  over 
Ro  We  note  that  x  is  an.  indeterminate  over  R0 


Definition.1:   Let  ye  E1   be  an  Indeterminate  over  R, 

1  n 


<2o2ol8> 


Then     u  =     Y  ®  Cn±  <f>  y1)      for     ui   €  R 


1=1 

4- 

n  €  Zq  is  said  to  be  a  polynomial  In  y  over  Ro  The 
set  of  all  such  expressions  is  denoted  by  R[y]„ 


From  the  definition  of  R[y]  and  the  previous 

0   1? 

observation  that   {x  9x  „x  „o00}   forms  a  basis  for  R1 


12 


we  can  conclude  that  R[x]  =  R.,  o  Also0  for  some  other 

indeterminate*   z,   over  R  we  can  see  that  R[z]  =  R[x] 

by  defining  the  function  X  i   R[z]  — — >  R[x]  where 
n  n 

u  =  V  ®  (u1  ©  z1)  £  R[z]  and  uX  =  Y  ®  (u±   $  xi)0   X  is 

1=0~  i=0 

easily  shown  to  be  an  isomorphism  from  R[z]  onto  R[x]* 


Definition?   Let  u  g  R[xL   Then  it  is  known  that  there 

<2o2cl9>     exists  an  n  such  that  for  all  1  >  np 

u.  =  0  where  u*   is  the  1-th  coefficient  of  uQ   If 

u  i  0C  u   is  said  to  have  degree  deg(u)  =  n=   If 

u  =  0   or  if  it  is  not  known  whether  u  =  0  or  not0   u 

is  said  to  have  virtual  degree  n.   The  degree  of 
0 

0  =  Y  (0  •  x  )   is  defined  to  be  =<*  , 

1=0 

Definitions   Let  n  =  deg(u)„  u  ^  0,   Then  u   is  called 

<2o2o20      the  leading  coefficient  of  u  and^Cu)  =  u  . 

n 


Definition?   Let  u  6  R[x]0   Then  u  is  said  to  be  monic 
<202o21>     lf/(u)  =  lo 

If  u0v  6  R[x]„   then  deg(u  ®  v)  <  max(deg(u)c 
deg(v))   and  deg(u  ©  v)  <  deg(u)  +  deg(v)o   Since  R[x] 
may  have  divisors  of  zero0  the  Inequality  in  the  latter  is 
necessary o  However „  we  will  be  mostly  concerned  with 
integral  domains,,  in  which  case  we  can  state 


13 
deg  (u  0  v)  -  deg(u)  +  deg(v) 

Not^t^on  Qonyention:   Up  to  this  point  the  vector  operator 
-^  symbols  have  been  kept  separate  from 

the  ring  operator  symbols  in  formulas.   The  normal  con- 
vention is  to  omit  or  replace  the  vector  and  scalar  opera- 
tors with  the  ring  operator  symbols  or  with  concatenation 
for  multiplication  when  no  confusion  results.   Since  we 
will  be  working  with  an  arbitrarily  large  number  of 
integral  domains  and  finite  fields  each  with  Its  own  set 
of  operators ,  we  must  adopt  some  of  these  conventions  for 
compactness  in  expressions „  but  since  each  of  the  finite 
fields  will  require  a  different  symbol  for  its  operators 
we  must  take  more  than  usual  care  in  making  these  notation 

conventions o 

K 

As  examples,  let  uj  ^  R[x3<   u  =   )  ©  (u,  0  x  ) 

M  k=0 


r 


(v     @  xm)c      Define     N  =  max(M0K) 


v  =      >   tt)     v 

L,  <=,       m 

N 

Then  u  ®  v  =     )  #   ( (u     ®  v   )   ©  xn) 

-.  L>  — .         n         n      ° 

n=0 

u  ©  v  =  y  ©  {(  y  $  (u.  ©  v j )  ©  xn) 

L„~~  Z^  km 

n=0         n-k-^m 

We  now  make  the  convention  that  the  scalar  product  operator 
symbol  ©  will  be  omitted  and  that  the  vector  operator  sym* 
bo3  p  will  be  replaced  by  the  ring  operator  O,   The  above 
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K  M 

formulas  now  reduce  to     u  =     \  ®   ^uvx   ) »    v  =     /  ®  ^vmx    ^ 

k=0  m=0 

N 
u  ®  v  =  V  ®  ( (un  ©  vn)xn) 
n=0 

K+M 

u©v=  y® (  y  ®(ukovm))xn) 

n=0   n=k+m 

In  addition  we  will  also  sometimes  omit  the  ©  symbol 
following  a  summation  symbol  when  the  summation  is  only 
a  formal  representation  of  a  polynomial „  For  example 
three  formal  representations  of  the  polynomial  u  would  be 

K  K 

U  s=  (Ug.UpUp.oo.oo)  =   \  ©  (ukX  )  =   \  U^X 

k=0  k^O 

However,  the  subterm   )  ®  (uv  ©  vm)   cannot  be  reduced 

/ 1      k.  m 

n=k+m 
to   V  (u^  ©  v  )   since  in  this  particular  case  the 

n=kVm 
symbol  V  ©  represents  an  actual  summation  using  the 

ring  operator  ©  0 


2o3  Polynomials  over  Integral  Domains,,  Finite  Fields 


Definition;  An  integral  domain  is  a  commutative  ring  with  a  unit 

<2o3,01>       ,   ... 

and  with  no  divisors  of  zero. 


Let  (1,+,°]   be  an  integral  domain.  Then  l[x]   is  the 
ring  of  polynomials  in  the  indeterminate  x  over  I. 
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Theorems   {I[x],+,«}   is  an  Integral  domain* 
<2o3o02> 

Definitions  A  field  is  an  integral  domain  in  which  every 
aJa    J  non-zero  element  has  a  unique  multiplicative 

inverse*   If  we  define  {F9©»0}  to  be  a  finite  field,  then 
by  the  above   {FCx],©,©}   is  an  integral  domain* 

Definitions   Let  p9a  €  F[xL  p   is  monic  and  deg(p)  >  1* 
<2*3°0^>     Then  by  the  division  algorithm  over  F[x]0 
there  exist  unique  q,r  €  F[x]   such  that  a  =  (q  ©  p)  ©  r 
where  deg(r)  <  deg(p)*   Then  define  a(mod  p)  -  r* 


Definitions   Let  p  €  F[x]s  p   is  monic  and  deg(p)  >  1* 
<2*3o05>     Define  F1  =  (a  |  a  €  F[x],  degCa)  <  deg(p)}* 
Two  binary  operations  on  F,   are  defined  as  follows 

1  1 

©  s  FXX  ?!     ' — ~>  ?2.°        a'b  €  Fl  B>     a  ®  b  =  a  ©  b 

1 

©  s  P1XP1  — >  Fx„    a  ©  b  =  la  0  b)(mod  p) 

1  1 
Theorems   {F199B®}   is  a  commutative  ring  with  identity* 

<2o3*06> 

Proofs   F^  c  F[xl*   F™   is  closed,  associative  and  commuta* 

-1       1 
tive  for  ®  and  O     and  contains  the  additive  inverse  of 

each  element  in  F.,*   The  distributive  property  follows 

from  the  fact  that  F,   is  closed  under  the  operations  © 


X 
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that  F[x]   is  distributive  with  respects  to  $  ,  0  and 
that   (a  ©  b)(mod  p)  =  a(mod  p)  ©  b(mod  p)   for  a„b  g  F 

Definitions   p  €  F[x]   is  said  to  be  irreducible  if  ^ 
<2c3o07?     a,,b  €  F[x]   such  that  p  =  a  0  b.   p  is 
said  to  be  a  prime  polynomial  (prime  element  or  prime) 
in  F[x]  if  p  is  a  monic  irreducible  polynomial  in 
F[x]  and  if  deg(p)  >   1. 

1  1 
Theorems   If  p   is  a  prime  polynomial,,   {F,6©0©}   is  a 

<2.3.08>  fleld# 

Proofs   Since  p  is  prime  there  exist  no  a0b  €  F,   such 
that  a  ©  b  =  p  =  0  (mod  p)  and  hence  F,   has  no 
divisors  of  zero.   Since  F^  is  also  commutative  and 
has  a  unit,  it  is  a  field 0 


Let  F  be  a  field  and  F-  c  F[x]  as  defined 
<2o3.09      in  <2o3005>0  Define  the  function  M  by  the 
following s 


M  s  F[x]  — ->  ?1     where  a  €  F[x]  ->  aM  =  a  (mod  p 


Th 


eorem°   M  is  a  ring  homomorphism  from  F[x]   onto  F, 


°Je  or  more  concisely  M  s  F[x]   onTo>  Fv 


11  2   2 

Theorems      Let      {R,©,®},    {B,,®,®}.    {R2»®»®}      be  rings »      Then 
<2o3o20>     ^    i    R     JL^     R^    T^    s    R^     aE_B>     R^     =>     T  = 

(TjTg)    i    R     -S->     R2 

Proofs   Let  O  represent  either  of  the  operators  ©  or  © 

1 
and  let  a„b  €  R«   Then   (a  Q  bjTjTg  =  (a^  O  bTx)T2  *= 

aTlT2  °  bTlT2  usinS  the  fact  that  T-,   and  T2  are 
homomorphisms  and  that  aT^bT,  €  Ri.  Hence  the  composite 
function  T  =  (T-.T-)   is  also  a  ring  homomorphism„ 

2*k     The  Modeled  Domains 

In  this  section  we  proceed  to  define  a  sequence 
of  polynomial  rings  over  an  integral  domain „   In  polyno- 
mial manipulation  on  a  computer  the  polynomials  with 
which  one  is  working  are  normally  in  theee  rings  and  there- 
fore we  call  them  the  "modeled  domains86* 


Definition g   The  ring's  I.  ,  1  >  0  are  defined  as  follows 
IQ  Is  the  ring  of  all  integers,,  an  integral 
domain  with  respect  to  the  operations   (+„°)o  For  1  >  0, 
let  I*  -i   oe  a  ring  with  respect  to  the  operations   C+„  v 
Define  1^  ■  Ii„1£x^]0   By  theorem  <2o2o10>  {Ilt)  +  o0}   is 
also  a  ring*   The  retention  of  the  operators  {-**0)      is 
justified  by  the  conventions  made  in  <2c2>0 


Lemmas   I*   Is  an  integral  domain 0 
<2.4.02> 


18 


Proofs   IQ  is  an  integral  domain  and  by  <2o3°02>  and  by  in= 
duction  on  i,   so  is  I.   for  0  <  i „ 


Lemmas  I±   c  Iy    0  <  i  <  j 
<2^<,03> 

Proofs   Strictly  speaking,  the  statement  is  not  true  since 
^1+1  ^s  a  vec^or  space  over  I.e   However,  an  embedding 

isomorphism  can  be  defined  which  maps  a  €  I .  into 

j 

(a90009  o  o  o )  =  ax.  ,  €  ^-«+-|o   The  normal  convention  is  to 
make  no  distinction  between  I.  and  its  embedding  in 
I.  ,„   a  convention  we  can  also  adopt  since  no  confusion 
will  result o  With  this  in  mind,  the  lemma  holds  trivially 
by  induction  on  h  =  j-i<, 


m 

Definitions   Let  a  6  I±    ?  a  =  V  afcxi»  am  +  °  for  °  <  1 
<2o4o04>                     k=0 

p.  oo  a  =  0 

0  0  *  j  -  i 

deg^(a)  =   V  m  0  <  j  =  i 

max(deg  .(a0)„ o o o 0  deg.(am))  0  <  j  <  i 

0  <  i  <  j 


Lemmas  a  €   I±     =>  deg.(a)  <  0,   0  <  i  <  j. 
<2o4o05> 

Proofs   I.  c  I.   and  by  the  discussion  in  <2o*K03>< 
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Definitions  Define  the  function  D^  ,  §  I±     —>     l^     as 


<2o^o06> 


follows s   Let  a  €  I« 8   if  I  >  00   let 


m 


a  =  I  akxi>  ak  €  X 
k~0 


1-1° 


aD 


»1 


fa 
m 


lk=0 


0  <  I  < 


<  1  < 


8  =1 

where  L  =  {b   b  =  aD,,  .  „  a  6  1*3 


ipj 


'i 


<2o^o0?> 

Proofs   Case  Is   0  <  i  <  js   by  <2o4.03> 


Case  2s   0  <  j  <  is   Let  h  =  i-jQ  Assume  that 

—    )     aipX4  Jul  o 


for  some  h  >  Q„  I,  cT,   Let  a  6  I 


1 


■i+r 


a 


m 


k~0 


-.1 


a,.D. 


ak  €  Ir   Then  ^±^^  L 

akDi,j  €  Xj  l 


€  I.      since 


dj  the  inductive  hypothesis  and      {I^o+o0} 

is  an  integral  domain  <>     For     h  =  10   a,DT    ,  =  akD=J    1   " 

=>1 

^-s  m   6  I^o  The  lemma  holds  by  induction  on 


Lemmas   D 
<2o^o08> 


H 


l.  J 


°>  I 


Proofs   Case  Is   0  <  i  <  js  D?  „  is  the  Identity  mapping, 

J-  0  J 

Case  2s   0  <  j  <  is  Assume  lemma  tea®  for 
-  K 

h  s  i  o  j  >  e0  Let  a0b  €  Ii+1?  a  -  £  akxi*l° 

k=0 
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b  -  I  V^i'  N  =  max(K*M)<   aDI+i„j  +  bDm,j  - 

N     m=0  N 

L(anDI>j  *  Vl^  -    I  (an  +  VD£j  "  (a  *  b)Dlil?J 
n=0  n=0 

=1 
since     D.     .      is  a  homomorphism  by  the  inductive  hypothesise 

K+M 

n=0  k+m=n 
K+M 

n=0  k+m=n 

=  (a  •  b)Diii,j 

For  h  =  0  the  lemma  is  true  by  case  1„   By  induction  on 
h  the  lemma  is  also  true  for  case  2» 

2<>5  The  Model  Domains 

Starting  with  a  finite  field  of  prime  order  we 
define  sequences  of  finite  fields  and  integral  domains  on 
adjunctions  of  the  prime  field «   Since  the  modeling  process 
will  involve  the  mapping  of  polynomials  in  the  modeled 
domains  into  these  finite  fields  and  integral  domains ,  we 
refer  to  them  as  the  "model  domains0 e   They  will  be  de~ 
noted  by  F.  and  J?  as  defined  in  the  following 
development o 

Definitions   FQ  =  {0019  2, » <> » .  t>P0~l}   where  pQ  is  a  prime 
<2o5o01>     integer  in  In„ 


Definitions   Let  @  and  ©  be  represented  by  Oo   Then 
<2o5o02>     define 

-1  -1 

O   §  IQ  X  ^o  —- ->  Iq  where  a„b  €  Xq  ->  a  O  b  -  a  O  b, 

-1 

More  specifically  a  @  b  =  a  +  bs 

=1 

a  ©  b  =  a  °  bo 

0  0 

O  §  FQX  F0  — >  FQ  where  a„b  €  IQ  »>  a  O  b  - 

(a  O  b)(mod  pQ)o 


0  Q 

Lemmas   {F0„@9®}   is  a  field, 

<2»5o03> 


Definitions  A  function  MQ  mapping  IQ  into  IQ  Is 

<2o5<>0*f>     defined  by  MQ'  s  IQ     —>     I    where 

aMQ  §=  a  (mod  pQ)   for  a  €  Iq  and 

I       -   {b    |    b  =  aMoo   a  €  Iq]     and     pQ     is  a  prime  integer « 

H 
Ti*£0£ej&s      M0    s    X0     onfo>     p0 

<2o5o05> 

Proofs  Given  in  most  texts  with  a  section  on  number  theory, 

Definitions  Define  E,,   to  be  a  finite  field  for  0  <  j 

°^°        with  respect  to  the  operations  ©  „  ©  o   Let 
p,  be  a  prime  polynomial  in  E.  ,[x.]o  Define 
d4  -r  deg.CpJ  to  be  the  degree  of  p„   in  the  lndeter~ 

mlnate  x°o 

5 
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Definitions  P.  =  {a  |  a  €  E,  ^x.],  deg,(a)  <  d,}   for  0  <  J 
<2o5o08> 

Definitions   The  operation  O0   representing  ®  or  ©  is 
<2o5<>09>     defined  as  follows  for  0  <  js 


O  s  P j  X  Fi  " — >     Yy   a°D  €  Fj   =>  a  O  b  =  (a  O  b) (mod  p 

Theorems   {F,,©,®}   is  a  field  for  0  <  j„ 
<2o5ol0> 

Proofs   By  theorem  <2o3o08>o 


Definitions  The  fields  E,  .,   in  definition  <2o5°07>  were 

0-Jo        unspecified  and  not  necessarily  related »  We 
now  define  E,   =  F,   for  0  <  jo   Consequently 
p.  €  F„  -iCx»]  and  the  operators  can  be  seen  to  be  consistent 

Theorem"   (Pixu^],©,©)   is  an  integral  domain* 
<2o5ol2> 

Proofs   By  theorem  <2o3o02>» 


Definitions  A  function  MJD  0  <  J,  associated  with  the 

°°  ^      prime  polynomial  p.  €  F.  ,[xj   is  defined  as 

follows  s   M,  s  F.  -.Cx,]  — — >  I    where  aM,  =  a  (mod  p.) 
j    j03  J.  j  p*  j  J 

for  a  €  F.   jlx,~}.      I   =  {b  |  b  =  aM.s  a  6  F..  -^Cx,]^ 

j 
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Theorem?   I   =  F^,  0  <  j. 
P  *  J 

<2.5-l4>    J 

Proofs   Let  a  €  F,  tCx,]-   By  definition  deg.CaCmod  pj)  <  d,( 

Therefore  deg.CaMJ  <  d,  and  aM„  €  F„Q  As  a  ranges 

over  P.  JxJ,,   all  of  I    is  mapped  outo  Hence 

I   cF„0  Let  b  €  P « o   Then     deg4(b)  <  d4o  Therefore 
P*  ""  3  j  J       J 

bM«  =  b(mod  pj  =  b  £  I    and  F«  el   » 

J  J         P  j        J  ""  P  j 

Theorems  M^  s  F^Cx^]   onTo~>  Fj<>  °  <  J 
<2o5ol5> 

Proofs  By  theorem  <2<>3<,10>o 

Lemma, s     P1  c  P,  ,[xJt  0  <  J 

<2o5ol?> 

Proofs  By  definition  of  E„  ,0  F,  and  <2o5oll>o 


Theorems 

F 

£ 

P. 

j  ° 

0 

<  i 

< 

3 

<2o5ol8> 

Proofs   Let 

h 

- 

j 

- 

io 

For 

h 

-  o0 

P*  =  Pj 

i    ,5 

c 

Fr 

Assume 

that  for 

some 

h 

> 

o, 

p 

'  ri 

c 

j 

» 

Then 

P    c 

j+1  £ 

F 

Ex 

rt1  b* 

<2»5ol7>  and  F^   is  embedded  in  both  F,[xJ+13  and  ] 

since  d.  >   1„  Hence  P,  c  F1+l  and  consequently  the 
theorem  holds  by  Induction  on  hQ 


£°rol3^rxs   Fa^  s  PjCxj^^  0  <  i  <  j 
<2o5ol9> 


Proofs   By  <2o5<>l8>  and  <2o5<,17>o 

Theorem?  F±   c  F,  c  P,  ^xJ  c  I.,  0  <  i  <  j 
<2o5o20> 

Proofs   The  first  two  relations  hold  by  <2o5°l8>  and  <2<>5ol9>< 

By  definition  FQ  c  I0„   implying  FqCx-^  -  Io^xl^  =  Il° 

Assume  that  F.  ,[x.]  c  I.   for  some   j  >  1.   Then 

F1  £  F1-l^xl-'  imPlies  F,  c  X,  and 

Fj[x„+1]  c  I<Cx1+l]  =  ^wj/   By  induction  on  j0 

P,  x[xJ  c  I,  for  0  <  jo 

Lemma s   a  €  Fi  =>  deg,(a)  <  0P  0  <  i  <  j 
<2o5o21> 

Proofs   Y  b  €  Fi,=i»  deg-(b)  <  0  since  the  embedding  of  P.  , 

into  F.  maps  all  elements  in  F,  .,   into  polynomials  in 

x.   in  F.  which  are  of  degree  0o  The  proof  is  completed 

by  <2o5o20>o 


Theorems   a  €  P,   =>  aM„  =  a 
<2o5<>22> 

Proofs   a  6  F.   =>  a  €  F,  1CxJ  by  <2o5<>20>  and 

deg,(a)  <  d.o   Therefore  aM,  =  a  (mod  p.)  =  aQ 

Theorems   a0b  g  P.   =>  aOb  =  aOb00<i<j 
<2o5o23> 

Proofs   Let  h  =  J  -  i„   For  h  =  0,   the  theorem  is  true 


trivially o  Assume  that  for  some  h  >  00  a  O  b  =  a  ©  b 
for  a„b  £  Fl0   a9b  <=  P,  c  FJ+1  by  <2o5°2Q>  and 

degj-fl(a)  <  00  degj+1(b)  <  0o   Therefore 


a   ^ 

b 

and     i 

leg 

holds 

o 

(a  6  b)  mod  p.,,  -   a  0  b  since  a  ©  b  € 
b)  <  Oo   By  induction  on  h„  the  le 


Proof?   By  <2o5°20>  and  <2o5°22>< 


The  sets  J?  are  defined  as  follows 
J 


<2o5o50> 


J  P1         tov     0  <   j  <  a 

!  J®.  ,[xj    for  0  <  a  <  j 


^JJ-1LXJ- 


Lemma  t      P  ..  c  J  ,  c  I  , 
<2o5o51> 

Proofs   Case  Is   0  <  j  <  a  by  <2o5o20> 

Case  2.8   0  <  a  <  . j  s  Let  h  •-  j  -  aQ  For  h  b  Q„ 
P,  -  J,  and  F.  c  J.  c  I.0   Assume  that  for  some  h  >  00 
F«  £  J«  £  I*«   L©t  a  €  F1+x»   Therefore 

a  €  F  [x.+1]  c  J'^Cx,  1]  s  J^,,  by  the  inductive  hypothesis., 
Hence  F^  c  J^.  J^  c  1^  =>  J*^  .  J^CxJ+1]  c  IjCz^] 
^1+1°  ^le  ^emma  holds  by*  induction  on  ho 


26 


Lemmas   Fi->l^xi-'  -  Jj9  °  -  a  <  J 
<2<>5°52> 

Proofs  F.,  c  J?  !     by  <2<,5o51>  and  hence 

ViCxj]  £  jj=iCxj]  =  J" 

Lemmas   J J  c  J*  0  <  1  <  J 

<2o5o53> 

Proofs   Case  Is   0  <  i  <  j  <  as   by  <2o5°20>„  <2o5ol8>o 

Case  2s   0  <  i  <  a  <  js   j£  =  F±   c  F^Ex^]  c  J^ 

by  <2o5°52>o 

Case  3s   0  <  a  <  i  <  js   Let  h  =  j  -  i o   For 
h  =  00  J?  c  J?o  Assume  that  for  some  h  >  00  Ji  c  J.  and 
a  €  J^o   Then  a  €  J  *  c  J^Cxj<fl]  =  J^r   The  case  holds 
by  induction  on  h„ 

Theorems   {J?,©9®}   is  an  integral  domain o 

<2o5o5^> 

Proofs   First  it  must  be  established  that   (®„$)   are  valid 

binary  operators  on  J.  X  J*» 

Case  Is   0  <  j  <  as   Ja  =  F,  and 

a„b  €  F.   =>  a  Q  b  =  a  O  b  by  <2o5o23>°   Consequently 

a  a  a 

{J  „©„©}   is  a  field  and  hence  an  integral  domain. 


Case  2s   0  <  a  <  j  s   Let  h  =  J  -  aQ   For  h  =  1 
£JaVle®9  J   is  an  integral  domain  by  <2o3»02>o   Since 
J?+1  =  J?txj^i^  fJa+h°®°  J   is  an  inteSral  domain  by 
<2o3<>02>  and  induction  on  hQ 


2o6  Further  Conventions  on  notation 

In  chapters  3  and  4  it  will  "be  convenient  to  have 
available  a  set  of  conventions  about  functional  notation 
which  will  help  to  reduce  the  size  of  expressions  and 
equations  involving  functions 0  We  therefore  define  the 
following  s 

Definitions   Let  A  be  a  set   (A)1  =  As  A11**1  «  (A)n  X  A„ 
««/  ss  I,  a   a  s  ififtoSi  oS«j  o  o  o  j  a  -»  y »  a^  c  a  J  o 

Definitions   The  component  extraction  operator  K.   is 
<2o6o02>     deflned  by  k   s  (A)n  — >  As 

aK^  s  a^  where  a  =  (a^..,,a^,,,,a  «)o 

Definitions   The  function  distribution  operator  is  defined 
<2o6o03>     by  q  .  (A)ny  F  _^  u  ^p)n  whepe 

P  a  {f  |  f  g  A  ~>  B}.  a  €  U)n,  f  €  F  =>   (a,f )e  s 

1&QX  0  Sl=j  I  p  o  o  o  O^L*,   -]  i  Jo 

Definitions   Let  f  s  A  — >  B9  a  €  Ao   Then  define 
<2o6o0^>     a(f)n  m   a  fop  n  £  0?  a(f)n  s  a(f}n-lf 

for  n  >  Oo 


Conventions   It  will  sometimes  be  convenient  to  treat  a 

<2o6o05> 

binary  operator  as  a  unary  operator  modified 
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by  the  second  in  the  pair  of  operands .   Let  f  s  A  ^  B  >  C, 

b 
a  €  A9  b  e  Bo  Define  af  =  (a,b)f „  The  following  are  some 

examples s 


(&Q9&19a2)   =  a  €  (A)3  =>  a(e)  K2  =  a£ffff 

f  o 
a(er  =  (a0ff0a1ffBa2ff) 


Definitions i   Ja  =  {j£  |  0  <  i) 

<2o6o06>      j  =  {Ja  |  „i  <  a}  (J-domains) 

I  =  {I.  |  0  <  i}  =  J=1       (modeled  domains) 
F  =  {Fi  |  0  <  i}  (F  domains) 

The  sequence   (p0»P15p?»o«o)   of  <2o5»0?>  will  sometimes  be 

referred  to  as  the  "prime  context  sequence" B 


Definitions   U  =  {Zq0~«>}0   The  symbol  «  has  the  following 

properties  i   Va  €  ZQ,  °°  +  a  =  °°;  °°  -  a  =  °°; 
a  <  °°o 
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CHAPTER  III   THE  RECURSIVE  MODELING  FUNCTIONS 

3ol  Introduction 

At  this  point  we  have  available  a  sequence  of 

modeled  domains   (Iq9I,„  0  <> » ,,1^  „  <> . « )  which  are  integral  do=> 

mains  with  respect  to  the  operations   (+0°)o  Based  on  a 

prime  integer  pQ  in  IQ  we  have  defined  a  sequence  of  F 

domains   (F00F10 „ „ » ,F., « o e )   each  of  which  is  a  finite  field 

generated  by  the  sequence  of  prime  elements   (PqoPto  °  <.  •  0p  *  „  o . 

where  p»   is  a  prime  polynomial  in  Fj..[x.3  for  0  <  j0 

The  ring  operations  on  the  sequence  of  F  domains  are  repre- 

0  0   11       j  j 
sented  by  the  operator  sequence   { (@„©) ,  (®„©)  „  o ° » „  (©, 


0  o  o  o  ,?  o 


Based  on  the  F  domains  sets  of  sequences  of  J  domains 
were  generated  one  of  which  is  represented  by 

(P0,F1M,,0Faol,Fa,Ja+lM,,0J.M,,)   where  FQ  through  Fa 


are  finite  fields  as  m  the  F  domains  and  J-0  5J„  .30°°° 


ra   ,a 
crt-l"  a+2' 


are  integral  domains  formed  by  successive  adjunctions  on  F  < 

The  corresponding  operator  sequence  is  then 

0  0   11       a-1  a=»l   a  a   a  a 

( ($„®)  o  (B0S)  o  o  o  o  „  (  e  o  G  )  ff  (e„6)  o  (9,®)', .  - . . ) 


The  object  of  this  chapter  is  to  define  sets  of 
functions  which  map  an  element  from  any  of  the  modeled  do- 
mains into  any  of  the  F  or  J  domains  and  which  in  addi- 
tion are  homomorphismso 
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3=2  The   C  Modeling  Functions 


Definitions  The  function  C    s  Ii 
<3«2o01> 


■>  I,  m      is  defined 


by  the  following?   I.  ,  =  {b   b  =  aC„  „  „  a  €  I,} 

Let  a  €  la  «   Then  if  i  >  0„   we  assume  a  has  degree  m 

m    x 

akxi  where  ak  ^  Ii„i •  Define  the  auxiliary  function 
k=0 
C.  .  §  I.   >  I.  .  where  I.  .  =  {c   c  =  aC.  . ,  a  €  I,} 

and  the  image  of  a  is 


aC 


i.J 


m 


I  (akCj-l«J-l)xJ 


k=0 
m 


^1/    ' 


L"®  (akci-i^} 


k=0 
aC 


i.j-1 


for  0  =  i  =  j 

for  0  <  i  =  J 

for  0  <  J  <  i 

for  0  <  i  <  j 


Then  define  aC.  .  =  aC.  .M. 


where  MQ  §  IQ  — — >  FQ;   aMQ  =  a (mod  pQ) 


M3  !  Pj-ltxJ]  — *  PJ!  aMJ  -  a(mod  V  °  *  3 


A  necessary  and  sufficient  condition  that  the  above 


definition  be  valid  is  that   I.  Q  c  IQ  and 

» 
*i  1  ^  Fi»i^xi-'  fo:r  °  K   «5  °  This  condition  is  shown  to  h 


old 
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in  Theorem  <3o3<>l6>o   In  addition  it  is  shown  that  1^^  ,  gF, 
and  that  C,  .   is  a  homomorphism  from  I,   into  F,o 

3o3  The  C  Modeling  Theorem  Proofs 

The  purpose  of  this  section  is  to  prove  the  fol= 
lowing  two  theorems s 

I±  Q   c  IQs  1^  ,  c  F,  x[xJ  for  0  <  J   <3<>3ol6> 


Cioj  °   Ij  "*"*  Pj   <3o3o25> 


The  above  two  results  are  the  only  ones  required  for  further 

developments o  If  the  reader  is  willing  to  accept  their 
validity  he  may  skip  the  rest  of  this  section  and  go  to 
section  3°**-° 


*      0  <  3   < 

J  0   J 

<3o3o01> 


Ijgmffi&s        a  <E   I,      ->     aC„     .   -  aC , 

J  i  0  J  ,  j 


Proofs  a  €  I*  ^>  a  €  Ii  ^  <2o^o03>o   Let  h  =  i  =  J, 
For  h  s  0„   the  lemma  holds  trivially*  Assume  that  for 

0  0 

some  h  >  09  aC,  ,  =  aC„  |0  a  €  I*   ->  a  €  I4,n   and 

i  o  j      Jo  J  0      l  l'*"-1' 

4e«1+1U>  =  0.  a<     =  I   J*  «aC;   )  -  «<   =  aC*   . 

k-0 

By  Induction  on  h  the  lemma  is  true» 


if 


<3o3o02> 
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t 
Proofs   CQ  0  is  the  identity  mapping  from  IQ  onto  IQ  and 

hence  IQ  Q  =  IQo   Let  a  €  IQo   Then  aCQ  Q  =  (aCQ  0)MQ  = 

aMQ  £  FQo   Therefore  ij  Q  c  FQ,   Let  b  €  FQ„   Then  b  €  IQ 

by  <2o5ol8>  and  bCQ  Q  =  aMQ  =  a  by  <2.5<>2^>o   Therefore 

OB 

Fo  -  Io0o° 


Lemmas    I    c  F^CxJ,  I~   c  F,   for  0  <  j 

<3o3o03> 

m 

Proofs   For  j  =  1„   let  a  6  I*  j  a  =  V  akx.,s  a,  €  IQo 

8      m  3  k=0 

aCl„l  =  I  (akC000)xl  6  F0[xl]  since  akC0„0  *  F0  by 
k=0 

8  8 

<3o3°°2>o  Hence  L  ,  c  F0£x-.  ]   since  I,  .,   is  mapped 
out  completely  as  a  ranges  over  I-0  Let  b  €  I-.  •  Then 

bCl  1  €  F0^xl^  and   ^bCl  1^M1  €  Fl°   But 

bC,  ,  s  (bC,  1)M,  6  I-,  -j   implying  that  1=,  ,  c  F,  „   Assume 

8  _    _  oo 

for  some  j  >  0  that  I*    .  c  F.  tLx^J  and  I.  .  c  F„0 

k 


Zk 
akxHl°    ak  ^   Il°      Then 

aCJ+l.J+l  -    I  (akc3,j)xj+l  €  Fjtxj+1]     since 
k=0 

oo 

aC.  ,  ^  I.  .  cF.  by  the  inductive  hypothesis „   Therefore 

JoJ         J  »  J   ""     J 

z'i*x,i*i  £  fjCxj+iL  (aCj+i„j+i)MJ+i  €  Fj+i  ■"* 
aCj+i.j+i  -  (aCj+i.J+i)Mj  €  ^l.j+r  Therefore 

oo 

1 1+1  .  -j  c  F„  ,„   By  induction  on  j   the  lemma  holds . 


Lemmas    I~^  g  Fj 
<3o3»04> 
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Proofs  A  combination  of  <3<,3°02>  when  j  =  0  and  <3o3o03> 
when  j  >   Oo 

Lemmas    c!  ,  =  C   ,   0  <  i  <  j 
<3o3o05> 

Proofs  Let  h  =  j  °  i  for  some  i  >   0„  Let  a  €  I,. 

Then  for  h  =  10   aC,  ,.-  =  aC.  .   and  consequently 

0 

C.  o  ,  =  C,  .   since  a  can  range  over  all  elements  in  I. 0 

o 
Assume  that  for  some  h  >  0„  C„  .  -  C*  . „  Then 

it j    it i 

0  c 

aCi  1+1  ~  aCi  1  ~  ^aCl  1^M1  ~  ^aCi  i^M1  by  the  in(iuctive 
hypothesis o   Since  aC.  .  €  F«   by  <3o3o04>  and  F»  c  F. 

by  <2»5ol8>  we  have  aC.  i  €  F„   and   (aC,  .  )M.  -  aCi  , 

by  <2o5o2^>o   Therefore  aC%    .+1  =  aCi  t  and  C1  .^  ^  (^  A 

By  induction  on  h9   the  lemma  holds 0 


Lemmas   b  e  F«  ^>  bC„  .  =  b 

j  j  t  j 

<3o3o06> 

Proofs   b  £  F  „  =>  b  €  I,  by  <2<>5o2Q>o   For  j  =  00 

b  6  FQ  and  bCQ  Q   =  bM(Q=  b  by  <2o5°24>0   Suppose  that  for 

some  j  >  00  b  €  F«   ^>  bC„  .  ■  bo   Then  let  a  6  F  s .-,  s 

m  J        Jo  J  j  "*"-•- 


a  s  iVj+i5  ffi  <  degj+i(pj+ih  ak  €  Fr   Then 

o  m  A 

aCj+l.J+l  s  kI0CakCj0j)xj+l  s  JQakx^l     by  the  inductlTe 

hypothesise      Therefore     aC.+1    „+1  €  F.Cx.^]     and 

o  m 

aCj+l„j+l  s  (aCj+lP j+l)Mj+l  s  I  akxj+l  ■  a  since 

ks=0 


m 
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<   deg,  ,(p.  J.   The  lemma  holds  by  induction  on  jo 


oo 


Lemmas   F^  c  1.^ 
<3»3<»07> 

Proofs   For  j  =  0,  FQ  =  I~  Q     by  <3o3»02>o  Assume  that  for 

oo  v     k 

some   J  >  0  FjCl^j.   Let  a  €  FJ+1;  a  =  £  akxJ+1? 

ak  €  Fr  m  <2e8j+i(pj+i)-   Then 

k=0 

m 


aCJ+l.j+l=  (  I(akcJ.3)xJ*l)MJ+l 

k=0 

m 

"  I  UkCJ,j)xJ+l  SlnCe  m<deS^l(p^l} 
k=0 


Therefore  FJ+1  c  I~+lpj+1°   The  lemma  holds  by  induction 
on  jo 

Lemma.?   F,  ^x.]  £  *<  i>   °  <  J 

<3»3o08> 

m 

Proofs   For  J  =  1,   let  a  €  FqCx^s  a  =  £  akxls  ak  €  F0 


k=0 


a  €  Ii   by  <2o5o20>  and 
m 


<i  =  I  <Vo.o»4  ■  IVi  =  a  €  :l.i  by  <J'3'06>' 

k=0         „  k=0 
Therefore  FqCx-^  Q  11  1"      SuPPose  that  for  some   J  -  " 

fj-iCxj]  £  xj0r   Let  a  €  fj[xj+i];  a  =  JLa***+lJ  &k  €  FJ 


k=0 


35 


Since     a  €   I,+1     by  <2„5.20>, 


aCj+l.J+l  -    I  (akCJ„J)xj+l  =    X/Vj+l  "  a     fey  <3«3.06>, 
k=0   ,  k=0  , 

Therefore     a  €  ^^t    «+x     a*"3-     Fi^xi+i-'  £  IjiJ.i    ^-n  °     The 
lemma  holds  by  induction  on     j 0 


_oo 


Ti'ftTlmft  °        ^  *    •*   —  F  jj 

J    9    J  J 

<3«3o09> 

Proof?  by  <3o3o0?>  and  <3o3oO^>< 


Lemmas    IJ0j  =  *J-l^x^f   °  <  j 
<3o3ol0> 

Proofs  by  <3o3o03>  and  <3o3o08>o 


o 
Lemmgs    Ji  j  =  Fjl°   °«i<^ 
<3o3oll> 

Proofs  Let  a  €  I1  for  some  I  >  0o   Then  aC,  .  s  aC,  „  €  F, 

by  <3o3<»05>  and  <3o3o04>o  Therefore  1^  ,  c  F10  Let  b  €  F^ 

Then  b  €  I-   and  bC„  ,  s=  b  by  <3<,3o06>. 

0  D  0 

bCj    .   =  bC.    .    -  b  6  Ii    ^     Therefore     F.    c  I.    . o 

o  _ 

<3o3<>12> 

Proofs   Let  h  =  I  =  j   for  some  fixed   j  >  0„   For  h  =  09 
^1  1  ~  Fi=.iLx^  by  <3o3olO>0   Suppose  that  for  some  h  >  0„ 

8  (1 

I±   j   =  F,   1LxJo      Then     b  6  1^     =>     bC±    ,   6  F  1=lCxj D o      Let 
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k 
keO 

r^  j^l    ' 

aCi+l  1s)®  ^aCi  1^  ^  F1-l^x1^  by  ^e  induc^ive  hypothesis 


a  €  I1+1s  a  =  £  ak*i+ls  ak  €  Ii°   Th« 


k=0 

0 

and  by  <2o5<>12>o   Therefore  Ii+1  «  £  ^i„i^xi-'0   ^et 

b  €  F1-l^xj^s  b  =  A  bkxV  ^k  €  Fj-1°  But  then  b  €  *j  -  Ii'fl 


by  <2o5o20>o   Therefore 

"kl.J  -  bcj„J  "  I  'Vj-lJ-l'^  =  I„Vj   =  b  *  <3°3°01> 

k=0  k=0 

and  <3<>3o06>o  Therefore  b  €  ^^.n  <  and  'i-v**^  -  Ii4-1  1 

implying  that  F,  ,[xj  c  I.  .,  «  a     The  lemma  holds  by  induction 

on  ho 


oo 


Lemmas    1^   =  F± ,  0  <  i  <  j 

<3»3ol3> 

Proofs   Let  h  =  J  -  i   for  some  i  >   0o 

Case  Is   h  =  Os   The  lemma  holds  by  <3<,3o04>. 


Case  2s   h  >  Os   Let  a  €  I*o   Then 


9 

aC 


i0j  =  aCi  i  €  Fi  by  <3°3o05>  and  <3o3»0*+>. 
aCi0J  =  ^aCi  J^MJ  =  ^aCi  i^M1  =  aCi  i  by  <2°5»24>0  Therefore 
aCl  J  €  Fi  and  Ii  1  -  Fi°   Let  b  €  Fi°  b  €  Ii  and 
bCl„i  =  b  by  <3°3o06>?  b  €  P.  by  <2o5<>l8>o 

bCloJ  =  (bCloJ)M  =  (bC1^1)M  =  bM  =  b  by  <2o5o24>„   There- 
fore b  €  I.  *     and  P.  c  I?  ,.   Hence  i"  .  =  P.   for  h  >  0<> 


Lgmmas   I10  =  IQ 
<3o3ol4> 
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Proofs   Iq90  =  J0  by  <3o3o02>o   Let  a  €  IQ  g  I.  . 

aCl00  =  aC01)0  "  a  €  Ii9o  by  <3o3o01>o  Therefore  IQ  c  I  Q  . 

Suppose  that  for  some  I  >  00  1^  Q   s  lQ0   Let  a  €  I 

aCi*l,0  =  L    ^CakCi90}  s  I   KCi  0}  €  X0  by  <2o5o02> 
k^O  k»o     "      u 

and  the  Inductive  hypothesise   By  induction  on  i0 


Lemmas   Iiffj  -  FJ9   0  <  j  <  i 
<3°3°15> 

Proof?    Case  Is   0  =  j  <  i .   F0r  i  =  0,  I 

<3o3o02>o  Let  a  €  1^  Then  ac"  0  €  IQ  by  <3o3ol^>  and 

o  ° 


'•  io()o  s  Fo 

D 


aCi00  =  (aCi90)M0  €  P0°   Hence  Ii9o  £  Fo° 

Let   b  €  F0o   Then  b  6  Ifl  c  I   by  <2o5ol8>  and  <2o4o03> 

bClo0  =  bC0o£}  =  b  by  <3o3o01>  and  therefore 


o 
bCi90  =  *bCi„0^M0  s  bM0  "  bo   Consequently  b  €  I?  n 

P  e-    I 
0  S  xi90° 

Case  2s  0  <  j  <  is  Let  h  .  ±        j  for  som 

j   >  Oo      For     h  =  0„    I 


JdJ   "  Fj     by  <3«3o09>o      Suppose  that 
for  soma     h  >  0.    ^   =  Fy     Let     a  €  PjS   a  =     |  ^k, 

ak  i  T>i}   m  <  degjipjh   a  6  PjC  I     s  It  S  I1+~  . 
<  i  a  m 

aCi^i.J  ■  aCJ9J  -  J^Vj-l.  j-i)^  -  j^  -  a  by 

<3o3»01>  and  <3»3o06>„   Therefore  a  €  F.   c  F  .  [x  1 

D  J       j°l    j 

aCi*l,j  ^  Pj-lt*j3-   Therefore  aC^^  =  CaC^^  =  aM 


J 


a 


by  <3o3o01>  and  <2o5<>24>0  Hence  a  €  l"   «  and  F  c  1°°' 
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oo 


By  Induction  on  h„  F.  c  I,  .  for  0  <  j  <  i.  Let  b  €  I,0 
bCi  i  e  Fi-l^xi-'  by  <3o3«12>  and  hence  bC±    .   =  (bC^  .)M  €  F 
Therefore  I?  *  c  F   for  0  <  j  <  io   Consequently  I.  .  =  F. 

J-  »  J       J  J-oJJ 

for  case  2» 


t  0 

theorems   (Ii  Q  =  IQ)   and  (1^    ,  c  F,-1Cxj]9   0  <  j ) 
<3o3»l6> 

Proofs    Case  Is   0  =  j  <  is  by  <3„3.1^> 

Case  2s   0  <  J  <  is  by  <3o3ol2> 

Case  3s   0  <  i  <  js  I±  .  =  F,  g  F.  c  F,  ^x,] 

by  <3o3oll>  and  <2o5o20>o 

oo 

Theorems   I„  .  c  F, 
1„J  ~  j 

<3o3ol7> 

Proofs    Case  Is   0  <  j  <  is  I*   =  F.  c  P.   by  <3o3ol5> 

1  3  J  J  J 


Case  2s   0  <  i  <  j s  I*   =  F±   c  F.  by  <3o3ol3>< 


Lemmas    Ci  i  =  ci  i   for  °  -  1  -  J 
<3»3ol8> 

Proofs   Let  h  =  j  -  i  for  some  fixed  i  >  0o 

Case  Is   h  =  Os   The  lemma  holds  trivially. 

Case  2s   h  >  Os   Let  a  €  L. 

aCi  j  =  ^aCi  j'M1  =  ^aCi  i^Mi°  But  aCi  i  e  Fi  by  <3«3»°9> 
and  F±   c  F.   since  h  =  j  -  i  >  0„  Hence   (aCi  ^M  =  aCi  1 

by  <2<>5o24>o   Therefore  aC.  .  =  aC  ,   and  C,  .  =  C   .  „ 
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Lemmas   b  €  I,  ^>  ^i*  =  bCj0i°  0  <  ^   <  1 
<3o3»19> 

Proofs   For  any  i  >  j  >  09  b  €  I* • 


bC„  .   =  CbC,  JM,  =  (bC,  ,)M.   by  <3<,3°l6>  and  <3o3o01> 

1  o  J         1« J    j         JoJ    J 

J  0  J 


<3o3°20> 

Proofs   Since  I°!  .  =  P.  by  <3o3o09>„  we  need  only  show  that 

J  0  J   „  J 


for  a0b  €  I.,  aC.  .  $  bC,  >   -   Ca  +  b)C.  . 

■*■      JoJ.     JoJ  Jo  J 

and  aC.  ,  0  bC„  .  -  (a  •  b)-C.  . 

Jo  J        JoJ  J  9  J 

For   j  =  0,  CQ  =  MQ  j  IQ  ^L>  FQ  by  <2o5°Q5><>   Suppose 

that  for  j  >  09  Q      .    §  I   -£->  F  .   Let  a  € 
K  3.J     J         J        M 


a  : 

r   k 

s  Z,  akxj+l" 
k=0 

ak  € 

I,  and  b  e 

IJ+1° 

b  = 

m 

mXj+lS 

bm 

€  I,.  Defi: 
3     N 
t-  b  ^  £  (an 
JS=0 

ne  N 

■  max(K9M) 0 

Then 

a  • 

+  V 

)xj+l>  a  °  b  s 

K+M 

-HI 

nsO  k+m- 

(ak 

n 

°  V): 

n 

xj+l 

aC 

A.- 
1+1 D 1*1   <* 

bCj*i, 

J+l 

K 

=  C  I  (akC1 

It         J+1 

.3,XJ+1)MJ+1  ®  ( 

] 

M 

)xm 

,HJ+1 

K  * 
=  (  Y  ©  (a. 

Cj0j)Xj+l)Mj+l   • 

c  £* 

m-0 

(bm 

m 

C    ) 

JoJ 

^tl)M 

3+1 « 

J  .         .  k     J 


k-0  ffi^O  ^  o  ->  °  Jo^ 


N  j         j 

n=0 
N  j 
=  (  Y©  (Can  +  bn^ci  i^xj+i^Mi+i  by  the  inductive  hypothesis 
n=0 

■  (a  +  b,Vi.w 


•cj*i.j+i  G  bCj*i.j*i 

K+M   j  j  . 

nsrO         k*m=n  "^^^  ana 


K*M 

n=0         k+ms=n 


=  (  £  4  C  £  ®  (Cak  -  bm]C  j))xj+1)Mj+1  by  the  inductive 


hypothesis 


K+M  j 

-  ( I  •  (t  I  K  •  V]cj,j,I5»i,,w  ^^  by  the 

n=0     k+m=n  inductive  hypothesis 

=  (a  °  b)Cj4l0^l 

H 
Therefore  C,+1  «+i  §  I1+1  "~ — >  F1+1°  By  induction  on  J, 

the  theorem  holds « 


CMaUaax'    ct^  :  iA   — >   pj9    o  <  i  <  j 

<3o3o21> 

Proof;  For  0  <  i  <  J,  C.  .  u   C.  ,  by  <303ol8>0  Therefore 

cicj  8  Ii   ~">  Fi  £  Fr 


kl 


<3»3o22> 

Proofs   For  i  -  00  CQ  Q  Is  the  identity  mapping  and  therefore 

is  a  homomorphisnio  For  i  >  0„   let  a„b  €  I*+i  s 
K  M 

a  =  X  akxi+l°  b  =  X  bkxi+l°  N  ~  max(K0M)?  akobk  €  1^- 

k=0  k=0 

o  g      I 

Assume  that  for  some  i  >  09  C.  Q  s  I.    >  X*  q 

aCl+l90  ®  bCi+l„0 

=  £®  (akCi,0}  ®  I®"(bmCi  0*   by  <2°5o02> 
k^O  m=0       ' 

N 


»     F®  Vxo^bnCi,0}      1V<2.*.02> 
ns=0 


N»l 
=     £®  ([an  +  \lc±   q)     *y  the   inductive  hypothesis 

n-0 
=  (a  +  b)C*+liQ 

C        =1     0  t  .8 

Similarly  aC1+1  Q  @  oCi+i  0  -  (a  °  b)C1+1  q0   By  induction 
on  i«   the  theorem  holds » 


<3o3°23> 

Proof?   G^  s  IjL  ^^>  F^1Cx^]  by  <3o3°12>o  Let  h  =  i  -  j 

for  some  j  >  0o  For  h  ■  00      let  a  €  I-}  a  =  Y  a^x^s 

M  k=0 

ak  €  Ijj]L  and  b  €  Iy   b  -  £  bmxm,  bm  €  Ij_r 

m=0 


hz 


Define     N  =  max(K.M) 

»         J-l        i 
aCJ5J     ®     bCJ,d 

=     I    ®     (akCj=l,j=l)xj     S       I    ®     (bmCj-lej»l)xj 
k=0  m=0 

-     I  ^    (anCj-l„j-l   ^  Vj-lJ-l'1"     *  <2-5-12> 
n=0 

N    .   -l 

"     I    •     (Can  +  bn]Vl,d~l)xJ     by<3-3o20> 
n=0 

i 
=   (a  +  b)C.    . 

J  9  J 

I  1-1  8  8 

Similarly  aC*  *  ©  bC,  .  =  (a  °  b)C.  ,„  Assume  that  for  some 

h  =  i  -  j  >  0,  ^  .  s  Ii  -*->  F  ^xJo   Let  a  6  Ii+1? 
K  M 

a  =   I  akxi+r  ak  e  Ji    and    b  €  xi+i»  b  =   I  Vi+r 

k=0  m=0 

b     <=   I,  0      Define     N  =  max(K„M) 
m         l 

8  ^"=1  8 

aCi+l0j         ®        bCi+l0j 

rH,  »       .     J-l      ?   H    ,  »       % 

=  J®  (akci„j)  •  Le  (Vi,j> 

k=0  m=0 

=     I    ®     (anCi0j     ®     b^   J      by  <3o3ol2>,    <2.5*12> 
n=0 

N  ,  -, 

=  )  ©  ([a  +  bv.^ci  4)   by  tne  inductive  hypothesis 
z_i_       n    n   1 0  j 
n=0 

=  (a  +  b)C°+loJ 


Similarly  ac"^^  ©  hC^^   -  (a  •  b)C1+lBj.   Therefore 
the  lemma  holds  by  Induction  on  ho 


Corollary t      C±^    s  I±     -B->  Fy      0  <  j  <  1 
<3o3o2^> 

Proof s    Case  Is   0  =  j  <  i 

Ci  0  s  Ji  ~>     V  M0  s  X0  "^  F0  by  <3°3o22>9  <2o5o05>o 

Hence  C1§j  -  C^  i  l±     -*-»  FQ  by  <2.3o20>o 

Case  2s   0  <  j  <  I 

<2o5ol5>o  Hence  C^  -  c"^  i  IJ>  J 


>  F,  by  <2»3-20>< 


<3o3o25> 

Proofs    Case  Is   0  <  I  <  j  s  by  <3„3o21>o 

Case  2s   0  <  j  <  i  s   by  <3o3<>2^>o 

3<,4  The  D  Modeling  Functions 

The  F-models  as  defined  by  the  C  modeling  fun©° 
tions  are  useful  in  themselves „  but  they  have  the  disadvantage 
of  destroying  the  degree  information  in  polynomials  <>   Co: 
sequently  we  define  another  function  from  the  modeled  doma, 
into  the  J  domains  in  such  a  way  that  the  degree  informa- 
tion of  selected  Indeterminates  Is  preserved,,  with  the  addi- 
tional requirement  that  the  functions  be  homomorphismso   It 
will  be  shown  that  this  is  accomplished  with  the  D  ffi'>i©liag 
functions o 


to 

Definitions   Define  the  function  D?1  .  :  I.   >  1^    ,     where 

J  i.j.oi  €  ZQS  ^  .  =  {b  I  b  =  aD1  ,.,  a  €  1^  as 

m 

akxl°  ak  ^  ^1-1° 
=0 


aC.  ,  for  0  <  i  <  a?  0  <  j  <  a 

1  0  J 

m 

k 


aDi.J 


I  (akDJ-l„3-l^J     for  0<a<i  =  j 

=0 

m  a 

][  ©   (akDi»l    i*  for     0<a<i?0<J<i 


k=0 

m 


k=0 
)a 


aD?   .   .  for     0<a<  j;    0<i<  j 


3o5  The  D  Modeling  Theorem  Proofs 

The  main  result  of  this  section  is  theorem  <3o5<»ll>« 
D,  .  s  I.   — ->  J.o  Again,,  the  proofs  presented  are  not  re= 
quired  for  the  developments  in  later  chapters  and  consequently 
can  be  skipped  if  the  reader  will  accept  the  validity  of  the 
above  theorem o 


00 

Lemma*   Iooo  =  Fo 
<3o5oOi> 

Proof 8   Iqo0  c  Fq   since  a  €  IQ     =>  aD^Q  .  aC0p0  €  FQ 

by  <3o3o02>o   Let  b  6  FQO  Then  bD^  Q   =  bCQ  Q   =  b  by  <3o3o06>< 

•'•  ro   £  10,0° 


Lemma  %        a  €  J*  =>  aD.  ,  =  a 
<3<>5o02> 


ks 


Proof 


Case  Is   0  <  j  <  a 


a  €  J?  c  I.   by  <2o5°53>  and  aD**  ,  =  aC.  .  =  a  by  <3o3c06>< 

Case  2  s   0  <  a  <  J 
Let  h  =  J  -  a  for  a  fixed  ac  For  h  =  0,  aD*J  .  =  a  by 

Jo  J 

case  lo   Suppose  that  for  some  h  >  0  that  b  £  J,      implies 

m  "~  » 

Q.  OL  V"1      k  CL 

bD.    ,   =  b0      Let     a  €  Jj+i?   a  =     )   akxHl?   ak  ^  J*°      TJlen 
a   €  Ij+1     and     a^  g  I    . 


m 


aD 


a 


j+l.j+l 


LkdLj)x 


m 


k=0 


j+1       JL^*^ 


k~0 


by  the  induct 
hypothesis 


=  a 


By  induction  on  h9   the  lemma  holds  for  this  case, 


Lejamas 
<3o5o03> 


Ta   Ta 


Proofs   Let  a  €  J 


J 


Then  a  6  I.  and  aD.  . 

J  Jo  J 


=  a  by  <3o5oG2>< 


Hence  a  €  I 


o  o   and  J ,  c  I   o 

J  9  J  J  J  0  J 


Lemmas   a  €  I-   =>  aD^  ,  =  &D*    .      0  <  j  <  i 

J  1»J  J  0  J  ' 

<3o5o04> 

Proofs    Case  Is   0  <  j  <  i  <  a 

a  6  I ,  c  Ii  =>  aD^  .  -  aCt  .  =  aC,  .  =  aD^  ,  by  <3o3°19>< 

Case  2s   0  <  j  <  a  <  i 
Let  h  =  i  =  a  for  some  fixed  a  and  j0  a€li  ->     a  6  I 
and  deg^a)  =  0  by  <2o4<,05>o   Assume  that  for  some  h  >  0 


=  i 
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,0i  T-vCl 


aD«  *   -   aD,  *<>   Since  a  €  *j*t»  degi+1(a)  =  0  and 

aDi*l  j  =  Z®  ^aDl  j^  =  aDi  J  =  aDJ  J   by  the  inductive 
hypothesis o  For  h=l0  i  =  a  +  1  and 

aD^,    b  aD^  ,  =  aD,  ,     by  case  10   Hence  the  lemma  holds 
for  this  case  by  induction  on  h„ 

Case  3°   0  <  a  <  j  <  i 
Let  h  b  i  -  j  for  some  fixed  j  „  For  h  =  0ff   the  lemma 


holds  trivially o  Assume  that  for  some  h  >  0 


aD,  ,  =aD.  ..   a  €  l°*+i      by  <2o4o03>  and  degi+-,(a)  b  0o 
Therefore  aD?.,  ,   =  aD?  .  =  aD?  .  by  the  inductive  hypothesis 

1+J-oJ       l  J       Jo  J 


Lemmas   a  €  J^  =>  aD^  =  a,   0  <  j  <  i 
<3o5o05> 

Proof  8    Case  Is   0  <  j  =  i  Directly  by  <3<,5°02>o 

Case  2  s   0  <  j  <  i 

a  €  J?  =>  a  €  I,  by  <2o5o53>  and  aD^"  ,  =  aD*?  -  ■  a  by 

<3o5o04>  and  <3o5°02>o 


Lemma  £   I^j  =  Jj 
<3o5o06> 

Proofs  J ^  c  I j  ,  by  <3o5o03>„ 

Case  Is   0  <  j  <  a 
Di  1  =  ci  1  =>  *i  1  =  XT  1  =  Fi  -  J<?   <3o3ol5>o 

X  0  J  J  S>  J  J  0  J  J  0  J  J  J 

Case  2s   0  <  a  <  j 

m 

Let  h  b  j  -  a.o   For  h  =  1,   let  a  €  Ia+1;  a  =  )   akxa+l 

k=0 


^7 


ak  €  xa  =>  aDa+l,a+l  =  I  (akD£,a)xa+l  €  JaCxa+l]  =  Ja+1 

k=0 

since  avD°J   €  J^  "by  case  lo  Assume  that  for  some  h  >  0 

m 

jj,j  £  Jr   Let  a  €  xj+i;  a  =  Iakxj+i°   Then 

m 

aDj+i,j+i  =  j^V^.j^i  «  JjCxJ+i]  =  Jj-i  sinoe  by 

the  inductive  hypothesis  a.D   .  £  J.0   By  induction  on  h„ 
the  case  holds  *     Hence  1^  ,  £  J?  for  0  <  a9    0  <  j 
implying  I*?  ,  =  J.0 


Lsnaa  *      d£<  =  dJ  1§    o  <  i  <  j 

<3o5°07> 

Proofs        Let     a   €  Ii 

Case  Is      0<i<j<a 
aD^    .   =  a^    .   =  a^   .    =  aD^  i     by  <3o3ol8>0 

Case  2  s      0  <  i  <  a  <  j 

Let  h  =  j  -  a  for  some  fixed  i  and  a0   For  h  =  1„ 

aD«  „.!  =  aD,  m  -   aD.  .   by  case  1»   Assume  that  for  some 
i » a+1     i  „  a     i  0  i   * 

h  >  00   aD.  .  =  aD„  .  »   Then  aD?"  .  ,  =  aD?  ,  =  aD?1  ,   and 
hence  the  case  holds  by  induction 0 

Case  3»   0  <  a  <  i  <  j 
Let  h  =  j  -  i   for  a  fixed  a  and  i„   For  h  =  0B      the 
lemma  holds  trivially,,   Assume  that  for  some  h  >  00 

aDi,j  =  aDi,i*   ^  aDl0j+l  =  aDi,j  "  aDi„i° 


Lemma .»    ij    =  j£„   0  <  i  <  j 
<3o5o08> 


48 
Proofs   D^j  =  Dj^  by  <3o5  =  07>  =>  l£  ,  =  lj9i  =  J*  by 


<3»5o06>. 


Lgmrn&s    1^  *  =  J*»   0  <  j  <  i 
<3<»5o09> 

Proofs   J.  g  I.  .  by  <3°5<>05>»  We  must  establish  that 
Ta     Ta 

Case  Is   0  <  j  <  i  <  a 
In  this  region  D?  ,  =  C.  .  <>  Hence  I?  »  =  I?  .  =  F.  =  J, 
by  <3<>3ol5>o 

Case  2s   0  <  j  <  a  <  i 

Let  h  =  I  -  a  for  some  fixed  a»  For  h  =  1„  let  a  €  I  +1  , 
m 


a  =  >  ©  (a.)  where  av  €  *a  1  =  Fi  by  case  1.  Therefore 

a  =  V  @  (ak)  €  F,  =  J^  by  <2o5o23>  and  <2o5o5^>o   Consequently 


k=0 

■  ^n  is  c  J.o  Assume  that  for  some  h  >  0,  I,  ,  c  J. 
a+1 o  J  —  J  i  o  j  ~  j 

m 


a  €  I? . -,  ,»   Then  a  =  Y  @  (a,)   where  a,  €  I?  .  c  J^  =  F. 

k=0 

m  •< 

by  the  inductive  hypothesise   Therefore  a  =  >  ©  (a,)  €  F.  =  J. 

k=0 
by  <2o5o23>,  <2o5o5^>o  Therefore  I?+1  *  £  J?  and  the  case 

holds  by  induction  on  h. 

Case  3s   0  <  a  <  j  <  i 

Let  h  =  i  -  J.  For  h  =  09  I*    ,  ■  J*     by  <3o5°04><,  Assume 

Jo  J        J 


a     *a   —      -a       Then 


that  for  some  h  >  0„  I.  ,   =  J.„   Let  a  £  *i+i  i° 
m  n 

©  (a,  )i   a,  6  I1  1  £  J.o   Therefore 


k=0 


2*9 

m 

a  =  I®  (ak)  6  Jj  by  <2°5°5^>  and  lJ+loJ  £  Jj°   Hence 
k=0 

I?  ,  S  J<   f°r  &11  cases. 

1  s  J        J 


Theorems   I.  ,   c  J. 
<3o5olO> 

Proofs     Case  Is   0  <  j  <  is   l£   =  J^  by  <3o5°09><> 

Case  2s   0  <  i  <  js   ij  .  =  jj  c  jj  by  <3o5*08>  and 

<2o5»53>o 


Theorems   d£   s  I±     -2->  J^ 
<3o5oll> 

Proofs    Case  Is   0  <  j  <  i  <  a 

Since  D.  .  =  C.  .   in  this  region  the  theorem  is  true  by 

1  9  J  1  9  J 

<3»3°25>o 

Case  2s   0  <  a  <  i  =  j 
Let  h  =  j  -  a  for  some  fixed  a0   Let  a9b  €  I 


K  ,  M 

k 


3+1* 

I  akxj*l'   ak   €lj     and     b  :      I  Vj+15    bm  €   V     Define 
k=0  m=0 

N  =  max{K9M)      and  assume  that     D?    ,    si.     -*— >     J?c 


aDJ+l  *  bDJ+l 

k=0  m-0 

N  N 

"   (  I*  (anDL>4n>  ®  (  X  +   (bnDL)xJ+l)by  <2c2o30><3o5ol0> 

n=0  n=0 
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N 


=  1}  (anDlj®Vj,j)xj+i     *<*-*-& 


n=0 

N 


=  T  ©  (Can  -i-  ^n)D?  Jx"^     by  the  inductive  hypothesis 


n=0 


la  ♦  b)D^ 

J   9    J 


K   ~  „     M 


=     L*(a^J.J)xJ+l®     I4   (bmDJ5J)xJ+l      *<2.2.30> 


k=0  m=0 

K+M 


=  £  ©  (  7  S  (akD^sJ  §  bmDjDj))xj+i  b^  <3o5ol0>0  <2„5o5^> 


n=0    k+m=n 
K+M  , 


=  £  ©  (  £  ®  ((ak  •  bn)Dj  j))xj+1  by  the  inductive 
n=0    k+m=n  '         hypothesis 


=  y  ©  ([  y  (a,,  •  b  )]D?  Jx1},,   by  the  inductive 
n=0     k+m=n  hypothesis 

K+M 


n=0  k+m=n 


=  (a  '  b)Dj+l,d+l 


Since  D^a  s  Ia  -JL->  j£  by  case  1. 

Da+l5a-H  %    Ja+1  JL>  Ja+lt   By  induction  on  h»   the  theorem 
holds  for  this  case0 

Case  3°°      0   <   j   <  a  <  i 

K  M 

Let     a,b  €   I1+1s   a  =     £  V^i    *k   €   ^     and     b  =     £  Vi+1! 

k=0  m=0 
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b  6  I* o   Define  N  =  max(K.M)   and  assume  that 
m    1 

Di  j  s  Xi  ~^~>  Jj   for  some  i  >  a* 


aDi+i,3§bDi+i,j 


vi*s£,>*l*<v?.,> 


k=0 
K+M 


=  Y    §  (  £  ®  (akDicj  0bmDi,j))  by  <3^5o09>9  <2.5o24> 


n=0    k+m=n 
K+M 


=  Y  ©  (  y  @  ((ak  °  "bm)D^  .))   by  the  inductive  hypothesis 
n=0    k+m=n 


K+M 


=  Y  ®  (C  Y  (ak  •  bm)]D^  .)   again  by  the  inductive 
n=0     k+m=n  hypothesis 


=  (a  '  b)Di+i,j 


By  a  similar  argument 


aD^   §  bD°"   =  (a  +  b)D?  . 


Let  h  =  i  -  cu   For  h  =  1,  D^  .  s  I±     — >  J^  since 
Da  „  s  I   — >  J   t»y  case  L   By  induction  on  h9 

Di  J  S  Ii  ™^">  Jj   f°r  °  -  J  -  a  *^  i* 
Case  4s   0  <  a  <  j  <  i 

Let  a9b  £  ^i+i   and  defined  as  in  case  3*   By  the  same 
reasoning  as  in  case  39  we  get 


52 


aD.  ,  ©  bD.  ,   -   (a  +  b)D.  . 


aD*  ,  ©  bD?  ,  =  (a  •  b)D* 

1  e  J        1»J  1«  J 


This  time  let  h  =  i  -  J.   For  h  =  1,  D*   ,  s  1^— >  J* 

since  D.  .  s  I,  >  J.  by  case  2»   Case  ^  then  holds  by 

j  J    J         J 

induction  on  h» 

Case  5?   0  <  i  <  j 

D*   =  D*   si   — >  J?  c  Ja  by  <3,5^0?>  and  <2  5  =  53> 


3^6  The  Recursive  Modeling  Function  and  its  Components 

Now  that  we  have  defined  all  the  domains  of  interesi 
and  have  proven  the  modeling  functions  to  be  homomorphismso  i" 
is  convenient  to  draw  all  the  results  together  by  looking  at 
them  in  a  different  manner  than  that  required  for  the  formal 
definitions  and  the  proofs  of  the  theorems 0   It  is  shown  that 
D,  .   can  be  expressed  independent  of   C„  ,  and  that  the  C 
functions  are  a  subclass  of  the  D  functions.   This  rela- 


tionship allows  us  to  call  the  D.  ,   function  in  its  new 
formulation  "the  recursive  modeling  function"  <>   Two  auxiliary 

functions  defined  as  component  functions  of  the  D,  .   func- 

i »  J 

tion  will  be  Introduced  for  the  purpose  of  describing  a  model- 
ing algorithm  when  the  polynomials  being  manipulated  are  in 
the  sum  of  products  form  rather  than  in  the  recursive  form, 


J. 


Lemma  t        C*    =  D^j,   0  <  j  <  i 

<3o6o01> 

=1         ! 
Proofs    Case  Is   0  =  i  ■  js   Let  a  €  IQ.  aDo9o  =  a  ~  aC090 

ft    k 
Case  2s   0  <  i  =  js   Let  a  €  I*,  a  =  )  akx,0 

ak  €   I^r 


j,j  -  I  («iP3:io-i,xj =  L(akcj-i.3-i)  =  aCis 

k=0  ksO 


aDJ"  V    (-    rJ" 


Case  3s      0  <  j   <  is      Let     h  =  i   -   j      and  assume 

1=1  ' 

that   for   some     h   >  0     Di    ,   =   C.    »o      Let     a   €   1jl+i» 

m 
a  =  kl0akXi+l9   ak  €   Ii° 


1-1  A   1=1  i_i  ft    1-1  o  » 

aDi+i.j  -  l  •  ifl.y  -  I  *  (Vi.j>  -  aCi*i 

k^O  k^O 


Therefore  the  lemma  holds  by  case  1  and  2  and  Induction 
on  ho 


Lemmas    Ci  j  =  Di  i»   0  <  i  <  j 
<3o6o02> 

Proofs  c[^   =  Cioi  =  D*fl  by  <3o3o05>< 

Lemmas    Di,j  ^  Di,i»   °  -  J  -  a 
<3»6o03> 


5^ 

Proofs  Case   Is      0   <   j    <  a,    0   <  i    <   j  s      D^        =   C  =  D^1    „ 

Case   2s      0<j<a0j<is      Let     h  =  a  -   j     and 

assume  that  for  some  h  >  05  D.  .  =  D;  . .   Let  a  €  I<Jl( 

"*     l » j    i  i,  j  1+1 


m 
a  =  L  akxi+l'  ak  €  Ii 


k=0 

aDiVj  "  I®  (akDi,J}  -  J®  <akDioj>  =  aDi+lej  *  <'2-5-23^ 
k=0  k=0 

and  the  Inductive  hypothesis <>   The  lemma  holds  by  induction 
on  ho 


Lemmas    C„  .  =  D^  7M , 

lo  3       lj   J 

<3o6o05> 

Proofs     Case  Is  0  <  j  <  is   D^~*M  =  cl     .M  =  C,  ,«, 

Case  2.  0<1  <  i,     Vi-)  .  J>i-\  .  ^   =  C^ 

DiIjMj  =  Ci9j  by  <3°5o0?>9  <3o6o03>  and  <3o6o02>o 


Corollary  s  D^        =  D^~*M, 
<3o6o06> 


Theorem;   The  modeling  functions  can  be  rewritten  in  the 

<3o6o0?>  .  ,,   . 

following  forms 

1  m 

Let  a  €  J°  ;   for   0  <  i   let  a  =  Y  a^*,  ak  €  J*1  . 

k=0 
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aD 


a 


i.J 


r* 


aDJ>3 


(1)      -1  <  a  <  i   <   j 


(2)      0   <  1   =   j    <  a 


S^-i.j-i^    (3)  °**<*- j 


k=0 
m 


I®   (akDi-l,j} 


k=0 


,aD 


i.i 


W      0   <   J    <  i4    -1   <  a 
(5)      0  <  i   <  j 


1  °° 


Proofs    Case  1  in  definition  <3o^o01>  (the  region  0  <  i  <  a, 
0  £  J  £  a)   can  De  replaced  by  the  following? 


aD 


a 


i.J 


aDL  -  Is  (akDi.i,J) 


k=0 

aD„  .  =  aD„  „ 
l »i      l 8 i 


»1  =  a,  0  =  i 
0  <  J  -  i  <  a 
0  <  j  <  i  <  a 
0  <  i  <  i  <  a 


by  the  use  of  lemmas  <3*6o01>  to  <3o6o06>o  After  the 
replacement  the  D.  .   functions  appear  as  follows s 

1  0  J 
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■X  =  a,  0  =  i  =  j 


aD^M  0  <  j  =  i  <  a 

J  j  J 


m 
aD" 


1 3  J  <   k=0 

aD^  i  0  <  i  <  j 


By  the  addition  of  the  a  =  =1  case  as  defined  in  <2o^o06> 
the  desired  formulation  is  obtained „ 


Definitions   The  modularization  modeling  function  D  °    is 
—  M« 

<3o6ol0>     ,  ..   ,     ^h>a  Th  . _-^   Ttt  *«,. 

J  defined  as  DM'  %  J.   - — —>     J,      for 

Mo  i         l 

~1  <  h  <  a  <  i„  0  <  lo    Let  a  €  J^o   a  €  J«    by  <2o5<>51>  so 
we  can  define  aD«°   =  aD.  ,  o 

Wj[         1  t,  i 


Definitions   The  substitution  modeling  function  Da.  .   is 

<3°6oll>     ,  „4        ^a      Ta Ta   _     _  _.  _ 

v  defined  as  D„.  .  s  J.   —— >  J„   for  -1  <  a0 

Ol  o  J       1  J 

0  <  j  <i„  a  <  jo   Let  a  €  J\  :   a  €  J?   so  we  can  again 
define  aDOJj  ,  =  aD.  .  „ 

Since  D,  .  s  J~.   — >  J   we  can  immediately  state 
the  following  lemmas? 

Lemma;    dJJ'01  s  J^  — >  J^„   =1  <  h  <  a  <  i  D  0  <  i  „ 
<3o6ol2>    i 
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Lemma  s    Dg     s  j£  — >  J*        -1  <  a,  0  <  J  <  i,  a  <  J , 
<3o6ol3>    ii;3 

*■»■     DS;a  Ds,    -  *5,  ,  D£;a>  -i  <  h  <  a  <  j  <  i. 

<3o6ol^>         *J      oJ   J 

The  purpose  and  Use  of  the  above  two  functions  can  be 

seen  more  clearly  if  we  consider  the  cases  in  which  a  -   h+1, 

m 
h      T*    k 
j  =  i-lo  For  a  €  J**  a  =  )  akxi  "the  modulariza'tion  and  sub- 

k=0 
stitution  functions  reduce  tos 


faMh+1  0  <  h  +  1  =  i 


aDM 
i 


,h,h+l     /  m 


I  (akDi!li^l)xi    0<h  +  l<i 
k=0 

m  h 
aD*      =    V  ©  (a,  )  -1  <  h  <  i9  0<io 

si0i-i    k£0    k 

The  modularization  function  states  that  xh+1  will  be  modular- 
ized and  that  the  degree  structure  of  x.   for  i  >  h+1  will 
remain  unchanged „   The  substitution  function  states  that  the 
coefficients  which  are  in  Ji=al  are  simply  summedo 

When  the  polynomials  are  represented  in  a  computer  by 
a  sum  of  products  form  successive  application  of  the  above 
modeling  functions  may  be  used  to  obtain  the  D.  .  model 0   By 
lemma  <3°601^>  the  modularization  and  substitution  functions  are 
commutative  which  implies  that  an  optimum  path  may  exist  which 
maps  a  modeled  polynomial  in  JT   into  its  model  in  J. o   This 
path  optimization  depends  in  general  on  the  prime  context 
sequence  and  on  the  structure  of  the  polynomial  being  modeled „ 
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3o7  Summary  of  Domain  and  Modeling  Relations 


a 


To  summarize  the  relations  among  the  I,,  J.  and 
F   sets  of  integral  domains 9    the  following  schematic  can 


be  studied 

0 

cc  . — >  »i 

0 

i 

2 

3 

k 

— -> 

i 

0    Jq1 

Jo 

h 

ss 

J-1 

0  <  i 

i      j;1 
i 

Jl 

ji 
Ji 

Jl 

= 

Fj 

0  <  j  <  a 

*    'i1 

4 

i 

4 

Ja 
J 

= 

JJ-l^j] 

-1  <  a  <  j 

3    ? 

^ 

*\ 

>\ 

>l 

k     Jlx 

>l 

*i 

4 

J3 

The  set  at  coordinate   (i0a)   can  be  thought  of  as  having  an 
adjunction  level  of  i  and  a  modularity  of  a,   where  the 
adjunction  level  is  the  number  of  indeterminates  in  the  domain 
and  the  modularity  specifies  the  indeterminate  which  is  the 
last  one  in  a  finite  fieldo   In  other  words „  the  integral 
domain  at   (i„a)   has   1   indeterminates  and  if  a  >  0„   the 
degrees  of  indeterminates  x,,x?ll.oo0x   has  been  reduced  by 
modularization  whereas  the  degrees  of  the  indeterminates 
x  -.pooooX.  are  still  possibly  intact  after  modeling., 
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The  above  schematic  provides  a  good  visualization 
of  neighboring  set  relations .   Define  the  following  symbols? 

C  -  Center  sets   Set  about  which  we  are  talking  -  (ego  (3^-)) 

L  -  Neighbor  set  to  the  left  of  C  -  (3,0) 

A  -  Set  above  C  -  (2,1) 

R  -  Set  to  the  right  of  C  -  (3,2) 

B  -  Set  below  C  -  (4,1) 

The  set  relations  then  are  given  by  the  following  % 

L   includes  C 

A   C  is  defined  by  the  adjunction  of  x,   to  A, 

C  includes  A  by  an  embedding e 
R   is  included  in  C 
B   is  an  adjunction  of  Co   B  includes  C  by  an  embedding 0 

Let  a  be  an  element  of  Co   Then  a  can  be  mapped  into  any 
of  the  k     neighboring  sets  by  functions  which  we  can  give 
the  following  names? 


c 

to 

L 

Reconstruct! on 

(a  €  C  => 

a  €  L) 

c 

to 

A 

Substitution 

Si,i-1 

c 

to 

R 

Modularization 

c 

to 

B 

Embedding 
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The  mappings  from  C  to  A  and  C  to  R  are  the  modeling 
functions  and  the  others  are  the  identity  function  for  the 
cases  we  have  so  far  considered  0   All  k-     mappings  are 
homomorphisms „ 

The  reconstruction  mapping  which  we  need  is  only 
the  identity  function  and  hence  information  lost  during  the 
modeling  is  not  recovered „   In  a  more  complicated  situation 
this  information  can  be  reconstructed,,   For  example,,  the 
Jl   to  JQ  modularization  is  commonly  called  residue  coding 

which  has  been  proposed  as  a  number  system  for  arithmetic 

(7) 
operations  in  a  computer  '    «      In  this  case  a  number  of  models 

(each  with  a  different  p~)  are  formed?  the  arithmetic  pro=> 

cesses  are  performed  in  the  model  domain  and  subsequently 

the  result  is  "reconstructed"  by  the  use  of  the  Chinese  re™ 

mainder  theorem . 

The  C  to  A  function  is  termed  a  'substitution" 

rather  than  a  summation  since  elements  in  Ji<=1   other  than 

the  unit  can  be  substituted  for  xi   to  produce  the  C  to 

A  function,,   If  we  substitute   3  €  J«  -<   for  x,   we  are 

essentially  finding  the  remainder  modulo  the  prime  polynomial 

(x.=3)   which  is  equivalent  to  the  modularization  from  C  to 

B«   But9  since  the  practicality  of  generalizing  this  to  a 

modularization  with  a  higher  degree  prime  polynomial  seems 

questionable ,  we  refrain  from  calling  the  C  to  A  sub- 

stitution  also  a  modularization,, 
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CHAPTER  IV  AUGMENTED  MODELING 

kd     Introduction 

The  D  modeling  functions  developed  in  Chapter  3 
have  the  ability  of  producing  models  in  which  the  degree  infor- 
mation of  selected  indeterminates  is  kept  intact 0   However,  it 
will  sometimes  happen  that  some  of  the  non-zero  coefficients  of 
a  modeled  polynomial  will  be  in  the  ideal   (Pq,P-,  *Pp»  <> «  <>  ,p  ) » 
the  kernel  of  the  D,  .  modeling  function.  When  this  polyno- 
mial  is  mapped  into  J.,  0  <  a  <  i  the  corresponding  coeffi- 
cients in  the  model  polynomial  will  be  zeroQ   If  the  leading 
coefficient  is  in  the  kernel ,  the  degree  of  the  model  polyno- 
mial is  lower  than  that  of  the  modeled  polynomial  a  circum- 
stance which  can  cause  an  algorithm  to  produce  erroneous  re- 
sults when  the  internal  branching  is  controlled  by  the  degree 
information  as,  for  example,,  in  the  division  and  g.c.d. 
algorithms, 

There  are  a  number  of  possible  solutions  to  this  pre- 
dicament „   The  first  and  most  obvious  is  to  use  a  very  large 
prime  for  pQ  and  high  degree  prime  polynomials  for  the  rest 
of  the  elements  in  the  prime  context  sequence.   This  reduces 
the  probability  of  producing  a  non-zero  element  of  the  D?"  . 

J-  9  J 

kernel  during  a  computation  to  an  insignificantly  small  value . 
It  may  then  be  assumed  that  the  degree  of  the  model  will 
always  be  the  same  as  the  degree  of  the  modeled  polynomial  and 
that  the  production  of  a  zero  in  the  model  implies  that  the 
corresponding  modeled  polynomial  is  also  zero.   In  some 


instances  this  approach  may  be  satisfactory ,  but  it  is  too 
inefficient „  often  dangerous  and  somewhat  inelegant  for 
general  use0 

A  second  solution  is  to  add  a  "virtual  degree" 
augmentation  to  the  model  in  order  to  be  able  to  detect  when 
the  leading  coefficient  of  the  modeled  polynomial  is  in  the 
kernel  of  the  modelings  the  condition  under  which  the  model 
will  be  called  "degree  reduced" «,   This  is  the  approach  which 
is  developed  in  what  follows «,   It  will  be  shown  that  the 
degree  information  can  be  treated  algebraically  in  a  manner 
closely  related  to  discrete  exponential  valuations    „   The 
virtue  of  this  approach  is  that  an  algorithm  in  the  augmented 
model  domains  is  capable  of  detecting  when  too  much  degree 
information  is  lost  to  continue o   Often  the  algorithm  can 
proceed  to  its  natural  termination  point  in  which  case  some 
information  about  the  result  of  applying  the  equivalent 
algorithm  to  the  modeled  polynomial  is  obtained,,   In  particular, 
the  "model  conclusive  g0c„do  theorem"  developed  in  this  chapter 
shows  that  the  goC0d0  algorithm  can  be  applied  in  the  model 
domains  and  zhat   the  algorithm  is  capable  of  determining 
conclusively  that  two  polynomials  are  relatively  prime  if  no 
degree  reduced  polynomial  remainder  is  encountered 0 

4o2  The  Augmented  Modeling  Domains  and  Functions 

In  this  section  we  define  a  set  of  augmented  modeling 
domains  which  consist  of  the  J  domains  to  which  a  "virtual 
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degree"  augmentation  has  been  added B   Once  these  domains  are 
available ,  sets  of  functions  are  defined  on  them  which  pre- 
scribe exactly  the  operations  to  be  performed  on  the  J  domain 
component  and  the  virtual  degree  component  of  an  element  in 
the  augmented  domain  when  the  function  is  applied^   Next  we 
introduce  the  concept  of  a  ' valuation  preserving  homomorphic1 
(or  VPH)  function  which  is  roughly  equivalent  to  stating  that 
a  function  having  this  property  commutes  with  the  augmented 
DM  modularization  modeling  function*   The  usefulness  of 
this  concept  stems  from  the  fact  that  composite  VPH  functions 
can  be  generated  by  argument  dependent  algorithms  composed 
of  VPH  functions  and  defined  only  on  the  augmented  model  do<= 
mains 0   In  other  words 0  an  algorithm  can  proceed  in  the  model 
domain  as  long  as  only  VPH  functions  are  met  and  at  each  point 
the  accumulated  composite  function  has  generated  polynomials 
in  the  model  domain  which  are  identical  to  the  model  polyno- 
mials which  would  be  obtained  by  applying  the  DM  modeling 
function  to  the  polynomials  in  the  modeled  domains  which  are 
generated  by  the  equivalent  accumulated  composite  function  on 
the  modeled  domains. 

These  concepts  are  not  developed  to  the  full 
generality  possible;  efforts  are  concentrated  on  obtaining  a 
set  of  definitions  and  relations  which  are  sufficient  to 
rigorously  prove  the  model  conclusive  g0c0do  theorem, 


6^ 


Definitions   (Extension  of  the  degree  function  cleg*  (ah 
<4„2.01>     a  €  ^  as  defined  in  <2o4,05>)   6^   s  j£  - — >  Ui 
a  €  J,   =>  a  €  I^o   Define 


deg^  (a)    0  <  i  =  j  <  a 
a6^   =  7  deg   (a)    (-1  <  a  <  i  =  j)   or   (-1  <  a,  0  <  j  < 
deg^  (a)    0  <  a?  0  <  i  <  j 


Note  that  deg*  (a)   was  defined  on  elements  in  I,   so  here 

we  first  map  the  elements  from  J.   into  I.   and  then  apply 

"l 
the  degree  function^,   Alsos  for  a  =  -1,  a6    =  deg,  (a) . 

U  U 

Definitions   Define  the  two  functions  ©  and  ©  by  the 
<4.2.02>     foilowlng! 


®  i    U  X  U  — >  UP  (asb)  ©  =  max(a0b)  =  a  ©  b 

u  u     u 

©  z    U  X  U  >  U,  (a,b)  ©=a©b^a  +  b 


where  U  was  defined  in  <2o6e07>. 


u  u 

Lemma t         {Uu©}   and   {U0©}   are  commutative  monoids  with 

,2o  3>  ag  the  additive   (@)   identity  and   0  as  the 

u 
multiplicative   (©)   identity* 

t(X        u.u.    ,   u   v  u  ,   u 
a,b0c  £  J±     =>  a  @  (b  ©  c)  =  (a  ©  b)  ©  (a  ©  c)„ 
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Proof i      Closure  and  associativity  follow  from  the  definitions 
and  since  max((a»b)sc)  =  max(a»  b,  c)  ■   =°°  is  the  additive 

identity  since  Va  €  J*,  a  ©  (~°o)  =  max(af=»)  =  a.   Similarly 

u 

a0O=a+O=a«   Distribution  follows  froms 


u    u    u 
(a  0  b)  ©  (a  0  c)  =  max(a  +  b„  a  +  c 

=  a  +  max ( b , c ) 

u    u 
=  a  O  (b  ©  c)  , 


Definitions   A  function   X  s  R  - — >  U  will  be  called  a 

<4  2  0k> 

valuation  homomorphism  from  the  ring   [Rt®}0] 

u 
into  U  if  va?b  €  R9   (a  ®  b)X  <  aX  ©  bX 

(a  Q  b)X  <  aX  O  bX 

Under  these  conditions  we  denote  the  function  X  by 

X  s  R  — >  U, 


Definition?   The  'Virtual  Degree*  or  V=modeling  function  is 

^'      ■>  defined  by  the  following? 

m 


Y 
k4o 


Let     a   €   J^    i   a  =     Y  a^xj 
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(1)  -1  «  a  <  i  «  j 

(2)  0  <  i  =  j  <  a 


aVj^  =  <   ae;1,  (3)   0  <  a  <  i  =  j 


WA  0  <  J  <  as  j  <  i 

I  ®  ^ak6i-l.  j}     WB  "1  -  a  <  j  <  i 

5)    0  <  as  0  <  i  <  j 

Define  V  =  {V^  ^  |  V^   »  Jj1  — >  Us  -1  <  a,  0  <  l,  0  <  j} 

_         _.oc     T»l   V    t» 
Lemma, '         V,  „  §  J\   — — >  u 

<4o2o05> 

Proof?   For  all  cases  but  (4)B   (-1  <  a  <  j  <  i)  s  aV?"  ,  =  a^T1,, 
— — —  —  i  „  j     i  5  k 

k  =  i0j      or     «>  = 

u 


(a  +  b)v£    .   =    (a  +  ^)S^Xk  <  maxCaS^^bd^)    =  aV^    .  ©  bV 


TIT  U 

a    •   b)vf    ,  =   (a   •   bid?1,    =  ae!1,    +  b6:Xv  =  avj1    .  ©  bV„a 

It  J  1»-K  19K  1,K  1 „  J  1 9 


J 


K  M 

Case    (4)Bs      -1  <  a  <  J    <  is      Let     a  =     T  akx^(,    b  =     V  bmxmf 

k=0  m=0 

N  =  max(K,M)0 

aVi,j®bVL=  Z^vlii.^5  Z®(vi-i,j) 

k=0  m=0 

n=0 
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n=0 

=  0 

a  +  b)v! 


> 

n=0 


k=0  m=0 

K+M  u  xi  u  -, 

=  I  •  <   I®  K6i-i;J  °Vl=i.j 

n=0         k+m=n 
K+M  u  u 

-  I     •   (      I®   ([ak    *    bm]6i-lJ}) 
n=0  k+m=n 

K±M  u  _  , 

>     )      9   (C      >      (a„        bra)]6.11     .) 
L>  Li         n  m        1=1  D  j 

n=0  k-t-m=n 

-  (a    «   b)V?    , 


Definition?   Let  b  €  J"1 .   Define  b  =  bV^"  ,  €  U   Let 

a  €  J*   and  define  a  =  (a.  a)   to  be  the  augmenta= 
tion  (or  augmented  model)  of  a*   a  is  called  the  D  com= 
ponent  (or  D-model)  of  a  and  a  the  V  component  (or 
V-model)  of  a   Note  that  a  =  aKQ;  a  =  aKx  by  <2,6  02> 
Define  j£  =  {&  |  a  =  (a,a)„  a  €  J^,  a  €  U} ,   Then 
±1   £  Jj  X  U, 


It  should  be  observed  that  the  V  component  is  not 
strictly  dependent  on  the  D  component.   An  augmented  modeling 
function  will  be  defined  which  assigns  a  value  to  the  V  com= 
ponent  of  every  a  €  «L  »   Then  sets  of  functions  on  L   will 
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be  defined  which  will  prescribe  the  exact  operations  to  be 

performed  on  :he  D  and  V  components  when  the  function  is 

applied-.   Hence  the  two  components  are  related  in  the  sense 

that  the  values  will  always  be  derived  from  the  same  set  of 

modeled  polynomials  and  by  the  application  of  a  parallel  set 

of  functions c   The  functions  defined  on  J„   will  be  such 

— 1 

that  a  and  a  will  always  be  related  by  the  inequality 

ct.     °  - 

aS«  <  £  a°   Because  of  this  property  a  can  be  referred  to 

1  9  1 

as  the  "virtual  degree"  of  a* 


Definitions   0  B  (00<=«>)0  1  =  (1,0),  ^1  =  (-ls0) 
<4»2<»07> 

Definitions   The  Augmented  modeling  function  is  defined  by 

<4o2o08>     ~a     _»1     .   Ta   _  .     ^  _-l   •  ^  TT    , 
D„  .  s  J.    — — >  J.e   Let  a  €  J4  t  a  €  U  and 

—1  0  J     —1  ~J  i 

a  -  (a,a)o   Then  aD,  «  =  (aD,  , ,  aV,  J:   Note  that  the  V 

^""1  0  J  1  D  J  1  «  J 

component  of  a  is  not  used  during  the  modeling-   Whenever 
a  modeling  is  performed  a  new  V  component  is  generated 
based  only  on  the  D  component  of  au 


Definitions   Augmented  modularizations   D^°a  s  J^  — >  J^9 
=1  <  h  <  a  <  1„   For  a  ^  £„  ^m   =  'a^M   t>a'° 


Definition;   The  augmented  addition  and  multiplication  func- 
tion @^  and  ©?  represented  by  O^  are  defined 
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as  follows?  Q£  s  (J^)2  — >  j£c   a,b  6  j£  =>  (a0b)o£  = 

~a        a    .  u  . , 
aQ^b=  (a  O  b,  a  O  b) 


Definitions   The  degree  function  on  J?  is  defined  by 
<^o2oll>     ,a     _a     .   TT   xa      -a 


*I  k  %   -i  ~"~>  -i°  For  a  =  (a„a)  €  j£«,   let 


Definitions   The  k8th  coefficient  function; 
<*K2°12>     _a   .  Ta 

m 
a  =  V  a,x1   for  1  >  0»   Then  the  augmentation  of  a.   is 

k=0 
a.  —  ^  a-.  „  a^  /  o 


-i.k  ~  A  ,_ 


09-<»)       (=1  =  a0  0  =  i)   or  (0  <  i  <  a) 
(aj^xj.c)     (=1  =  a9  0  <  i)   or   (0  <  a  <  i) 


where     c  =  0     if     a,     >  Ox    c  =  =«     If     a,    =  ~eo< 


Definitions   Leading  coefficient  functions 

<^o2ol3>     ao.   „  Ta  ___^  _a   T  .     _  ,a  .    v 

£*    o   J*   — — >  J„0   Let  a  €  J^o  k  =  aK,  =  a 


f  (09-*»)      if  k  =  -oo 
a^i  =  ■{    a 

Definitions   Reduction  functions 

<4o2ol4>      a   /Ta*2     .   Ta   T  .     .  r    Ta 
P<  s  (J*  )   — >  J  0   Let  aDb  6  «L 
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a 


b„         n  (c,a=l)   if   0  <  b  <  a9  b  =  b6    0 

ajg,  =  (a0b)p,  =  <  o 

-i        All  other  casess   (b  <  0),  (a  <  b) 

t«j    <  b. 


where  c  =  [((b*f)  0?  a)  ©?  ((-1)  0?  (a*?)  @?  b  g°  x^)]K, 

— —  1    ~ 1  "~   —1  —1    — *=-  1    — 1       1    1     V 

and  h  -  a  °  "6o 


Definitions   Residue  modulo  bs   H?  s  (J?)2  — >  J?o   Let 


<^o2<»15>     „  K  g.    Ta 
1,9  b  €  «L  < 


-i    v-i'        ~i 


b 
(a5b)R^  -  &(p£)  ,     h  =  a  -  b  +  1 


Note  that  the  function  is  well  defined  even  for  h  <  0  by 
<2o6o04>o  b  is  called  the  divisor  of  the  residue  function,, 

Definitions  f±  m   {f£  |  f^  s  (J^)n  — >  (J^)m0  -1  <  a  <  i} 

f,   is  called  a  valuation  preserving  homomorphic 

(VPH)  set  of  functions  on  J?  and  J*}     if  for  a  6  (J?)n 
and  -1  <  h  <  a  <  is 

nh  D a         nh , a 
1        a     h   1 

Define     VPH^c0L  =    {£      |  X±      is  a  VPH   set   of  functions   on 


-i 


[?     for     -1  <  h  <  a  ^   i} 
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Lemma?     f,g  €  VPH^' a  =>  fg  <E  VPH^a 
<4o2ol7> 

Proofs   Since  D^f'a  s  J^  — >  j£  by  <3*6*12>  the  D  compos 

ent  equivalence  follows  immediately <>  The  V  component 

equality  follows  since  aD» f  &..  =  a* 

Lemmas    ©r  0±   €  VPH^ a 
<^,2d8> 

Proofs      Let     a„b  €   J^6        and     O     stand  for    ©     or     ©  0 

(^^.bD^Q?  -    ((aD^'aObD^a)5    (a  O  b) ) 
i  i  i  i 

h  h   a        o  u    . 

=   £(a  O  b)D;Jt,(\    (a  Ob))      by<3  =  6„12> 

i 

=    (aQ*  b)D^a 


Lemmas  %     k   g  VPH^t0L,       -1  <  h  <  a  <  i 

<4o2»19> 

m  h 
Proofs      Let     a  €   J*;    a  =    (a.ah    a  =     Y®   (akxi^ 

k=0 


Then 


m 


a 


Y©   (a,D!j  a    *-k 


ak  =    (ak,ak),    aD    »a  =     £©   ^M^  )xi 


a^     Dhncx  _    ,        0      )Dh.a  _    f(     Dh0a   }    0      } 

aDha  a       _    (   Dh5a  &)    a       _    ((      Dh,a   .    0      } 
aiiMi  ii„k  "    laL,Mi    taJ£isk  ~    UakDMi^1Jxi^c; 


where     c  =   0      if     a,    >  0     and     c  =  »°°     if     a,    =  °°° 
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Lemmas  ^   €  VPH*"a 

<^.2o20> 

h  c  h  cu  * 

Proofs      Let     &  €   J**    k  =  a  =  aK,  .      aD^"    K,    =  a  =  k= 

ha  ^ 

k  =  =>*>     ->     aDJJ9      =  a  =  0o      For     k  =  -«\ 

a/hDh'a  =   0l£'a  =   0 

^m'^1   -    <aD£>a.-o>    =   0  j 

i  "x  i 

■o  i    ^  n  h_h  ,  a  h     _,h  ,  a  _h ,  a.  a  ^h  ,  a  a 

For    k  >  o,  a^  D       =  a*lik2„      =  aEM    2C1|k  =  a2M:  *r 

i  l        i  i  i 

b  ~"Mi 

<4,2o21>  11 

Proofs      The   condition     0   <  b&^       =  b  =  bD{j°a6  guarantees   that 

b^*  £  0     and     t^sCt^  ^  0:      Let     a  €   jj* 

b  ^M± 

Case  Is   a  <  b  ^>  a(pj)  -  a,  fl£JJtCX(pJ)   =  aD^a 

Case  2?   a  >  b.   The  'c9  function  in  the  £^  definition  <4,2.-l^> 
is  composed  of  functions  in  VPH„ '   and  hence  is  also  in 
VPHh:Cl,   The  fact  that  the  leading  coefficients  in  b  and 
bDM    are  not  zero  guarantees  that  the  coefficient  of  xa  in 

1*1-4  1 

c  will  be  zero  and  hence  that  the  virtual  degree  can  be 

reduced  m  both  domains, 

vh.  a 


— Mi 


eorems       (ifjDJJ'  a  =  D*»a(R*)        for     b   €  J £ 


JJti 
<4„2<>22> 
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0  <  b6*  ,  =  b  =  bD^tCX6?  .  .  ,  n 

b  m\ 

Proof?   Under  the  given  conditions   (p£)DJJ,a  =  D{j,a(P^)   0   Since 

IL   is  a  composite  consisting  of  only  £.   functions,  the 

theorem  holds,, 


4,3  The  Model  Conclusive  g.cd.  Theorem 

We  have  available  now  a  sufficient  number  of  functions 
defined  on  the  augmented  domain  to  generate  polynomial  remainder 
sequences  in  the  model  domains  and  to  develop  the  model  con= 
elusive  goCd,  theorem  with  which  we  are  able  to  show  that  two 
polynomials  are  relatively  prime  with  respect  to  a  specified 
indeterminate  by  performing  operations  only  in  the  model  domains 
A  number  of  examples  are  given  to  demonstrate  the  power  of  this 
theorem  and  to  illustrate  some  problems  which  might  be  en° 
countered  when  applying  the  theorem  in  practical  cases 


Definitions  Let  &„bDr  €  jj,  -1  <  a  <  I,  0  <  i.  Then  (a0btr! 

is  called  a  polynomial  remainder  triplet   (p.rjt.i 

a  °     "    a 

over  J.   if  a6,  >  b6  >  r&  and  there  exist  b,cL„  a_     g  J.  „ 

Jj  =  (b  90)0  a  =  (a  ,0)   and  b   6   =  a  6  -  0   such  that 

b  0  a  =  (b  ©  q)  ©  (a  ©  r)  „   It  follows  that  a  =  b  ©  q. 

Note  that   6  Is  used  to  represent   6„  „ 

1 , 1 


L   =  '.EqiLioomL.)   is  called  a  polynomial  remain* 
der  sequence  (pc.rso)  if  each  triplet 


7k 

'£m<£nn.r£nrt.2)  is  a  p°rttc  over  ^i  for  °-m-n=  2t  £ 
is  called  a  normal  p.rcS,  if  r_ ,  n  6  =  r6  -  1  for  1  <  m  <  n< 

r  is  called  a  complete  p0rcSc  if  lJ^-q   =  ° 

The  above  definitions  are  consistent  with  those  in 
references  (5)  and  (6)«   It  is  to  be  noted  that  the  V  com= 
ponents  were  not  used  in  the  definitions  and  hence  a  p,r,t, 
and  poToS,  could  have  been  defined  only  on  J.  • 


Theorems   (Euclidean  Algorithm)  If  r  is  a  complete  p0rcSc 

<4  3  03> 

aJ"    J        then  r  ,   is  an  associate  of  the  g0c„dc  of  rntr,o 
n^x  u  x. 

a 
Proof?   Since  J.   is  an  integral  domain,   r  =  0  implies 

r  t   divides  associates  of  r   for  0  <  m  <  n  -  1;   Hence 

n-1  m        —  ~ 

the  goCodo  of  rn,r,   must  be  an  associate  of  r  ,  «> 

0  1  n-X 


Let  t-   (lmJ,,,...!.)   and  r_  ,  n  -    ( r  ,  r_.  . ,  )  R , 


a 

K)5±-lso  v^n}      """  ^m+2  "  'Tn^+l'-i' 

J  Then  r  is  a  p^r.s.  if  and  only  if  rQ  >  r-,   and 


r  =  r  6  for  1  <  m  <  n.   We  call  r  the  p(r,st  generated 
mm         —  — 

by   (r()lErR1)( 


Proofs   Under  the  condition  r«  >  r,   and  r  =  r  6, 

0  —     1  mm 

r_   =    (r^(,r,  )R„;    r0  =  rn    -   1  by  definition  of  the     R4      function 

~*c  ~"U    ~"X    "~1  <c  X  ~ "1 

and      (r^ „r,  „r0)      is  a  p:r,tc  The  condition     r     =  r   6     guar- 
— >U  —i  —2  mm 

antees   that     r~  =    (r^rOfL  *  By  induction  on     ni;,      we   see   that 
^£o°~l°~l^      generates     r0 
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If  r  Is  generated  by   (E^r^IL),  rQ  >  r^  by  the 

definition  of  a  generator  of  a  pcT.s.  and  by  the  definition  of 
Re   If  for  some  k»   ^— v-?'— k-1'— k^   is  a  P°r»to  and  rv  ~  ° 
or  rk&  <  rk£,  (£k.!-£k)Sj  =  £fc_i'   Hence   (xk.i«£k»Ek+i)   is 

not  a  p»r.to  Therefore  if  r  =  (rQ,£,,...,rL)   is  a  p<,r„So 

o 

and      (rA, o o « f r    er    .t )      is  not,    the  condition     r     -  r   6     holds 
—0  —n  ~-n+l  m         m 

for     1  <  m  <  n. 


Definitions   r  =  (rQ9r,  „  «, » .  ,r  )   is  called  a  conclusive  ptrts 
<^o3o05>     if  r^  _  0  and/or  r^  _  0c  By  the  prececLing 

lemma  and  the  definitions  of  R. ,   a  conclusive  p.r.s.  must 
also  be  normal c 


The  above  definitions  and  lemmas  hold  for  polynomials 
in  the  modeled  and  model  domains  since  no  restriction  was 
placed  on  the  value  of  arj   the  modularization  level.   The 
important  problem  to  be  solved  is  to  find  a  relation  between 
the  two  which  allows  us  to  make  a  conclusive  decision  about 
the  modeled  polynomials  while  performing  operations  only  on 

their  models  in  the  model  domain.   The  specific  question  we 

=1 
wish  to  answer  is  whether  two  polynomials  In  J.   are  re= 

i 

latively  prime  with  respect  to  x.s 

Theorems   Let  En*— 1  €  J7  p  rQ  >  r-,  =  r,6  >  0:   Then 

v^o'-l'-i  '   generates  the  p0r0So   r  =  (r^r,,, „ <, „ 0rN) 


where     N  >  2.      Let      s     =  rD~1,a  €  J?e      Then     s0   >  L      and  if 

sl  ~  sl^i   1°    ^n9  — 1°—  l^     generates  the  p<,r0s  = 
s  -   (sA9St p o o o 9s    )     where     n  <  NQ 

b 

Proofs   (S1)  €  VPH^°a  if  0  <  b  =  b6h  .  =  bD^°a6?  ,   by  <^02o22> 

=~~ ~" ~""       ""1  1  °~" ""  1  ,  1      Tl.   "~1  o  1 

Hence 


«s0>Sl)H^=  <i0l£1'a.r1l£L*'X 


i 


iro^R-1^ 


"  ^M, 


=  s2 


l9a 


s06?  -  <  r^"1,   by  <4.2&17>,   By  induction  on  m»  (snv  sn  ,R?) 
generates   s  and  s.,  =  s  =1  for  m  <  n.   s„6,  .  =  0  or 

Q.       "• 

s  6,  „  <  s   t =1  and  consequently  the  p.rrS  terminates . 
-=n  i,  1    n=l  H     *     ^ 

s  =  r  D~     implies   s6.  ,  <  r  6~  ,  =  r<   Hence  n  <  N  = 
— n   — n^M,      ^      — n  i  ..  i  —  — n  i ,  i  — 

Theorems   (Model  Conclusive  g „ c  d ,  Theorem) „   Let  rss  be  the 
<4o3o075  p$7iSo  of  the  previ0US  theorem  <k*3o06>   and  let 

g  =  gcC.do  (rQ1)r-j)o   If  s  is  conclusive 

1)  with  s  =  0,  s  ^  0  then  g6°   =  0c 

2)  with  sn  >  0„  sn  =  0  then  g6~  1   <   s   -. 

Proofs   Case  Is   sR  =  0„  sn  ^  0s   By  theorem  <^3.06>  n  <  N 

-1       =1 
and  hence  n  =  N  and  rN  =  0o   Therefore  g6„  .  ■  rw6i  i  = 

~1  "1  a 

since  rKT6.  ,  =  -«>  would  imply  r  =  0  and  s  =  r  DM  4  '  =  0 
jm  1 0 1  n  n    n—n. 


which  contradicts  the  case  hypothesise 

Case  2s   s  >  0P  s  -  0%      r«  »,   is  an  associate  of 

-1  =1      •=! 

the  goCodo  of  rQ  and  r,   in  J.  o   Hence  r^^i^i  ±  ~   s^i  ± 

and  rN  =  Oo  But  n  <  N  by  theorem  <^o3o06>  implying 

sn-l  =  rn4  -  rN=l  =  rN=-l6i9.i° 

Stated  less  succinctly  the  model  conclusive  g0CocL theorem  states 
that  if  a  model  p0r0s«  is  conclusive  and  the  last  remainder  is 
not   0„   the  original  modeled  polynomials  rQ  and  r.,   are 
relatively  prime  in  the  indeterminate  x.o   The  second  condi= 
tion  states  that  if  the  last  term  in  the  p„r<>So  is  0„   the 
goCodo  of  rQ  and  r^  must  have  a  lower  degree  in  x   than 
the  degree  of  the  second  last  term* 


The  model  conclusive  gcCodc  theorem  has  much  prac- 
tical significance  since  it  will  show  that  two  polynomials 
are  relatively  prime  when  the  por.So  is  conclusive  much  more 
rapidly  than  methods  in  which  a  g0c0d0  algorithm  is  applied 
directly  to  the  modeled  polynomials 0   To  get  a  better  grasp 
of  the  implications  of  this  theorem  we  can  consider  the 
following  set  of  non-trivial  examples? 


ty orations   The  polynomials  will  be  represented  in  sequence 

<j^o3   10> 

form  as   explained  in  appendix  A, 


78 

Examples   Find  the  GCD  of  rQ  =  1  +  x?c,  r,  =  3  +  x-^ 

°^°      •*       g  =  goCodo  (r09r,)o   To  illustrate  a  conclusive 
PoToSo  in  which  s  =  09   we  use  p«  =  5° 

rQ  =  1   0   1   ,2 

rx  =  3   1      ,1 

1   2   0   ,1 
z3  =  0         ,0 

r~  =  0  implies  that  g6~  _  <  10   We  repeat  the  calculation 
using  pQ  =  7° 

rQ  =  1   0   1   02 

£l  =  3   1      »1 
14      ,1 

r2  =  3         ,0 

This  shows  conclusively  that  the  two  polynomials  are  rela= 

2 

tively  prime o   The  problem  for  pQ  =  5  was  that   1  +  x, 

factors  into   (3  +  x-,)(2  +  x1)   modulo  5„ 


Examples   The  following  polynomials  which  appeared  in  reference 

Wo  page  588  illustrate  a  number  of  problems  which 
may  be  encountered  and  at  the  same  time  show  how  effective 
the  model  conclusive  g,c8d,  theorem  is, 

a  =  (~x3)y6  +  (5x2)y5  +  (x3+x2~7x)y4  +  (~2x2~2x+3)y3 
+   (x3~2x2°2x*l)y2  +  (x2+6x+5)y  ♦  (~x3~3x2~3x~D 
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b  =  iVa 

Let  x1  =  x8  x?  =  yD   Then  the  polynomials  are  represented  by 


a  = 


* 

0 

1 

2 

3 

0 

..l 

-3 

-3 

-1 

1 

5 

6 

1 

0 

2 

1 

-2 

=2 

i 

3 

3 

-2 

-2 

0 

k 

0 

-7 

1 

1 

5 

0 

0 

5 

0 

6 

0 

0 

0 

=1 

b  = 


0 

5 

6 

1 

0 

1 

2 

=4 

-4 

2 

2 

9 

=6 

=6 

0 

3 

0 

-28 

4 

4 

k 

0 

0 

25 

0 

5 

0 

0 

0 

-6 

Let  the  prime  context  sequence  be   (pQ  =  5t>Pi 
Then 


X-n  —J.  9000/0 


aS 


2.2 


* 

0 

1 

0 

4 

2 

1 

0 

1 

2 

1 

3 

3 

3 

3 

4 

0 

3 

5 

0 

0 

6 

0 

0 

3 
1 
0 
0 


0 
1 
0 
1 
0 
4 


3 

2 

3 

2 

3 
2 
3 


<2 


ft 

0 

2 

0 

0 

1 

1 

0 

I 

2 

1 

1 

2 

2 

4 

4 

4 

0 

3 

0 

2 

4 

k 

h 

0 

0 

0 

0 

5 

0 

0 

0 

4 

£0 


I 

2,2 

*1  =  ^252 


((202,3.^0,0,^),6) 


((2„l,>2oOp094)05) 

2 
£2  s  (£o'£l}£  =  ^®  ((2,002,2,0)04) 


SinCe  M2.2 


3  ^  r?  -  4„   the  result  is  inconclusive „   Now 


we  use  the  prime  context  sequence   (pQ  =  110P-.  =  x,  +  10 


$    o  o  o  s   o 
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aD*^  =  ((3B1B901006„5,10)06) 
~202  =  ^1"7,81>2P395)95) 


^0  =  -aD2.2 

= 

8 

10 

2 

1 

5 

6 

1  06 

Xi  =  9©  (b£2j2) 

= 

9 

8 

6 

7 

5 

1 

o5 

(r^r^P 

= 

8 

1 

5 

6 

9 

1 

>5 

0 

£2  =  ^£oo-l^-°-l^- 

2   o 
r2  =  3  ©  r2 

= 

10 
8 

4 

1 

10 
8 

10 
8 

4 
1 

,4 

ter^ 

= 

9 

0 

5 

10 

8 

,4 

£3  =  ((£1»£2)£o£2^- 

= 

0 

3 

7 

1 

o3 

(£29-3^ 

= 

8 

1 

5 

1 

o3 

0 
£4  =  ( (£2»£3)£o£3)£ 

£4  =  5  ©  r^ 

= 

8 
7 

9 

1 

9 

1 

.2 

.2 

(£3»£i|)£ 

= 

0 

7 

6 

,2 

£5  =  ((£3»£i|.)£»£^)£ 

s 

2 

1 

.1 

(£5"£i|)£ 

= 

7 

10 

•  1 

£*  =  ((£«ror,l)p9r,,)p 

B 

9 

,0 

The  result  is  conclusive  and  indicates  that  the  g0c0d0  is  inde< 
pendent  of  yQ 

If  we  try  to  use  the  D  ,   model  to  check  if  the 
goCodo  is  independent  of  x-.  0   an  inconclusive  result  will 

always  occur  no  matter  what  the  value  of  pQO   This  is  caused 

3 
by  the  symmetry  of  the  polynomials  in  the  coefficients  of  x£ 

which  will  generate  a  zero  coefficient  during  summation,, 

£i  ■  p^2  1  wlH  De  degree  reduced  and  consequently  cannot  be 

used  as  a  divisor  in  the  reduction  function,. 
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If  It  is  still  desired  to  continue  in  the  model  do- 
main r  two  alternatives  remain.-   The  x-,   and  x?  indeter- 
minates  can  be  interchanged  by  a  substitution  operation  and  a 

prime  sequence*   (Pq^Pt  =  x,  «  p,o„o)   where   3^1  can  be 

1  2 

used  In  the  modelings?   rQ  =  &D  and  r,  =  bD2  _,   If 

this  still  does  not  destroy  the  symmetry  which  causes  the 

cancellation,   P-i's   of  higher  degree  in  x-,   can  be  used: 

Alternatively,,  since  it  has  been  established  that  the  g.ccL 

is  independent  of  y0   the  pair-wise  goCod*0s  of  the  y 

coefficients  can  be  calculated e   This  discussion  should  pro° 

vide  sufficient  evidence  to  Indicate  that  there  exist 

significant  problems  in  determining  the  best  strategy  to  be 

used  after  an  inconclusive  p.r.s-  is  encountered  in  the  model 

domain o 

Example t      As  a  final  example,  we  consider  a  36  term,  2  vari= 

<!(,  3  12> 

J  able,  total  degree  7  pair  of  polynomials  whose 

coefficients  are  decimal  digits  extracted  from  a  table  of 

random  numbers,   Collins  (6)  demonstrated  that  of  the  h   goC.d 

methods  he  considered  only  the  reduced  p.r«se  algorithm  was 

capable  of  showing  similar  pairs  relatively  prime  In  a 

reasonable  amount  of  computation  time.   The  prime  context 

sequence  is   (pQ  s=  ?9P-i  ■  x-,  -  !•.••)• 
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a  = 


* 

0 

1 

£ 

1 

& 

i 

6 

2 

Os 

8 

4 

6 

7 

8 

8 

6 

1 

Is 

8 

8 

4 

2 

8 

9 

2 

2s 

8 

6 

0 

2 

3 

5 

3- 

6 

7 

1 

5 

0 

4s 

5 

5 

5 

0 
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9 

3 
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7? 
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2 
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£0  =  ^2,2 

£o  =  4  ©  4 
i 
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.1 
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The  result  is  conclusive  and  Tq  £  0     implies  that  the  g.Cod, 
does  not  contain  a  term  with  x_o 

It  should  be  noticed  that  if  the  two  polynomials 
do  have  a  g.cd,  of  degree  1  or  greater  the  algorithm  in 
the  model  domain  will  always  be  either  inconclusive  or  have 
a  final  remainder  of  zero  in  the  sequencer   s   1   will  be  a 
model  of  the  actual  goC0d*  but  having  available  only  the 
developments  derived  previously  in  this  chapter  we  are  not 
able  to  calculate  the  actual  gcCodo  from  this  model „ 

In  a  polynomial  manipulator  the  model  g.cdo 
algorithm  should  be  used  as  a  "filter"  which  is  able  to 
catch  most  polynomial  pairs  which  are  relatively  prime 0 
depending  on  specific  values  of  elements  in  the  prime  con- 
text sequence o  The  polynomial  pairs  which  do  survive  this 
filtering  will  very  likely  have  a  non-trivial  g0c0da  At 
this  point  the  filtering  can  either  be  repeated  with  a  much 
more  refined  prime  context  sequence  (e0go  pQ  is  larger)  or 
goCodo  algorithms  such  as  discussed  in  reference  (6)  can  be 
applied  in  the  modeled  domain. 

kok     On  the  Use  of  Augmented  Modeling  in  Elimination 

Let  &■.     €  I   and  let  the  equation  a,  =  0  be 
represented  by  e,o   The  system  of  equations 

E  =  Cek  I  ek  °  ak  =  °   °  —  k  —  n  ~  "*"}   then  has  n  equations 
in  m  unknowns  which  are  in  general  not  linearly  related c 
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The  elimination  method  is  the  orderly  application  of  the  re- 
duction operator o   P»   to  the  equations  in  the  system  so  as 
to  attempt  to  reduce  the  problem  to  that  of  finding  the  zeros 
of  an  equation  b  =  0  where  b  g  L.   For  the  case  in  which 
E  is  a  linear  system  of  equations  the  method  reduces  to  the 
commonly  used  Gaussian  elimination  in  which  the  system  E  is 
represented  in  matrix  form. 

The  operations  used  during  elimination  are  similar 
to  those  used  in  finding  the  g0c0d<,  of  two  polynomials  and 
consequently  we  can  expect  to  be  able  to  generate  some  con- 
clusive results  by  operations  performed  only  in  the  model 
domains  *      As  an  example  we  will  demonstrate  the  use  of 
augmented  modeling  as  a  method  for  checking  the  consistency 
of  the  system  E  when  E  is  linear  and  overspecif ied  in  the 
sense  that  n  >  m<» 

LiQir     e  =  16a(  G-i  06nj  6«/  j  e-^.  s  a-^.  =  u v   a-,  c  ^-o 


eQs  x~  +  2x2  +  3X-L  +  4  =  0 

e-^  2x~  +  4x2  +  5X-L  +  6  =  0 

e2?  x~  +  9x  +  6x1  +2=0 

e_s  x~  +  3x3  +  5^  +  1  =  0 


e: 


The  above  is  a  system  of  4  equations  in  3  unknowns  in  which 
the  initial  ordering  was  made  arbitrarily „   Since  we  will 
have  to  be  able  to  Interchange  the  ordering,,  we  define  the 
following  operator: 
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Definitions   Sn  .  s  (A)n  >  (A)n;  S*  ,      induces  a  permuta- 

1  4  J  -L  0  J 

<4C4„01>     tlon  Qn   (A)n  guch  that  fQr 

a  =  (a0,a,,,o.»a  J  6  (A)  ,  aK„  =  aS„  .K   and  aK.  =  aS^  .K^ 
where  the  component  function  K  is  defined  in  <2o6o02>« 


For  convenience  we  make  the  following  notation  conventions 


e  ^  «2  —  >6n(6«i  i  e  /->  ^  —   n  l  '  9     ~   o  o5°3^  ~  ^SrvoSiu^^vSo 


The  system  of  equations  are  represented  in  the  matrix  form 


1  2  3  4 

2  4  5  6 
19  6  2 
13  5  1 


First  we  generate  the  D>   .  model  of  the  system  with  pQ  =  7» 

The  V  component  of  the  augmented  models  appears  after  each 

comma;  only  the  V  components  which  will  be  used  in  the  elimin- 
ation  are  shown. 


1,1  2,1  3,1  4,0 
2  -  pD°   .     2«1  ^X  5,1  6,0 

1,1  3,1  5,1  1.0 


The  application  of  the  elimination  method  on  the  system  results 
in  the  following  set  of  system  states 1 
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0 


1,1  2,1  3,1  fc.O 

t     -    (2  o    roUc   0,0  0,1  6,1  5<0  degree  reduced  in  x. 

e   -  ve^^  ^f  <   o,0  0,1  3,1  5,0     "       "       "  n* 

0.0  1,1  2,1  4,0 

,   -,  1.1  2,1  3,1  M 

I  -  (e   e    S3   )•   °>°  1'1  2>1  4*° 


2 
*1 


0S0  0,1  6,1  5,0 
1,1   2,1   3,1   *K0 


2  _  ^2        ,At\.  0,0   1,1   2,1  4,0 

0,0  0,0  3,1  5,0 

o„o  o,o  6,1  5„o 


-~  (%),iv~253(-^c  0,0  o„o  3,1  5,0 


i              1,1  2,1  3,1  4*0 

e  -  (2     1  (0))-    °»°  lsl  2sl  ^° 

e-  iSo|2tS^l£JJ-    o„0  050  3,1  5,0 

0e0  0,0  0,0  2,0 


4 
The  last  equation  in  e  states  that  2  =  0(mod  7)i  a  contra- 

1  11 

diction »   e.  is  degree  reduced  in  its  equations  e-,      and  e?* 

Since   e-,   should  be  used  as  the  divisor  in  the  reduction, 

the  method  would  fail  since  £  would  not  be  in  VPH.   By 

performing  an  interchange  of  equations  with  the  operator 

3    2    1 
S«  p(,  e-,=  e~   is  used  as  the  divisor  in  the  reduction  which 

3 
generates   et   As  a  result  the  entire  sequence  of  operations 

performed  on  e  are  all  in  VPH  and  consequently  theorem 

<4-2  17>  applies.   Since  equation  e~   is  a  contradiction ,  a 

corresponding  contradiction  holds  when  the  equivalent  sequence 

of  operations  in  the  modeled  domain  is  evaluated.   In  other 

words ,  the  above  contradiction  in  the  model  domain  implies 

that  in  the  modeled  domain  the  rank  of  the  augmented  matrix 

of  the  linear  system  E  is  greater  than  the  rank  of  the 

coefficient  matrix. 
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It  is  of  interest  to  note  that  the  sequence  of  opera- 
tions applied  in  the  model  system  creates  an  "operation  trail" 
based  only  on  decisions  derived  from  the  model  system   When 
this  same  trail  of  operations  is  applied  to  the  modeled  system 
it  is  guaranteed  that  the  divisor  of  any  reduction  operation 
will  not  be  zero  and  consequently  that  the  reduction  will  be 
performed  effectively o   This  in  turn  implies  that  as  long  as 
the  model  operations  are  valuation  preserving  homomorphic 
functions,,  all  sequencing  decisions  in  the  elimination 
algorithm  can  be  made  independent  of  the  information  in  the 
modeled  system,, 

Though  a  further  study  of  the  applications  of 
augmented  modeling  has  not  been  made  it  is  expected  that  its 

use  can  be  extended  to  nonlinear  systems  and  to  the  reduction 

(2) 

of  systems  of  differential  equations  as  discussed  by  Brahns , 

Also0  under  the  VPH  proviso  the  solutions  generated  in  the 
model  domain  are  models  of  the  solutions  in  the  modeled  domain 
and  consequently  if  a  number  of  model  solutions  using  different 
primes  are  generated,,  a  reconstruction  based  on  the  Chinese 
remainder  theorem  which  generates  the  solution  should  be 
capable  of  being  defined  for  at  least  some  classes  of  linear 
systems o   In  reference  (9)  where  Moses  discusses  the  prac= 
ticality  of  solving  systems  of  polynomial  equations  by  elimina- 
tion it  is  brought  out  that  the  main  problem  is  the  explosive 
growth  of  coefficients o   Since  augmented  modelings  can  be 
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defined  so  that  the  coefficients  are  elements  of  finite  fields, 
this  growth  is  completely  under  control.  However,  the  problem 
that  now  presents  itself  involves  determining  whether  or  not 
sufficient  information  is  still  retained  in  the  model  to  draw 
conclusive  results  from  it* 
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CHAPTER  V  THE  PRIME  CONTEXT  SEQUENCE 

5ol  Introduction 

In  previous  chapters  the  existence  of  prime  context 
sequences   (Po»PicPp»°t">)   was  tacitly  assumed-   Since  the 
primes  in  these  sequences  determine  the  properties  of  the 
modeling  we  will  briefly  look  into  their  characteristics  and 
generation,,   The  prime  integer  pQ  plays  the  most  Important 
role  since  the  rest  of  the  polynomials  in  the  sequence  are 
defined  on  adjunctions  for  which  pQ  determines  the  prime 
subfieldo 

5>2     PO 

The  prime  integer  pQ  determines  the  number  of 

elements  in  the  prime  subfield  and  therefore  can  be  used  to 

control  the  amount  of  storage  required  for  the  models a      To 

maximize  the  information  content  in  relation  to  the  number 

of  bits  of  storage  in  a  binary  machine „   pQ  should  be  of 

n 
the  form  2     =  m  where  n  is  the  number  of  bits  and  m 

is  the  smallest  odd  number  which  optimizes  some  other  aspect 

of  operations  in  FQo      These  other  aspects  cannot  be  specified 

independent  of  the  particular  method  and  equipment  used  to 

perform  the  operations „   For  example „  if  the  operations  are 

performed  by  subroutines  in  a  computer  with  an  unalterable 

structure,  not  much  of  an  increase  in  speed  can  be  achieved 

by  choosing  different  values  of  m   However,  if  the  hardware 
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in  the  computer  is  alterable  in  such  a  way  that  it  can  for 
example  perform  ones  complement  13  bit  arithmetic „  addition 
and  multiplication  modulo  8191  can  be  performed  very  rapidly 
whereas  a  choice  of  p0  =  2   -  31  would  require  much  more 
time  to  perform  the  same  operations  given  the  identical 
hardware  configuration e   This  then  is  an  example  of  hardware 
being  naturally  tuned  to  a  particular  pQ» 

If  the  circuitry  in  the  computer8 s  central  pro- 
cessor is  not  rigidly  specified,,  it  may  be  possible  to 
change  it  dynamically  so  that  it  will  always  be  tuned  to  the 
particular  pQ  being  used.   This  in  turn  implies  that  addi- 
tion and  multiplication  modulo  pQ  can  be  performed  in  about 
the  same  amount  of  time  required  to  perform  integer  addition 
and  multiplication  respectively »   That  this  is  possible 
should  be  apparent  by  observing  that  the  end  around  carry 
in  ones  complement  arithmetic  corresponds  to  the  reduction 
of  a  number  modulo  2n  -  1.   If  the  prime  is   2n  -  m9   the 
end  around  carry  digit  is  replaced  by  the  addition  of  m  to 
the  least  n  bit  positions  of  the  number ,   Showing  that 
this  results  in  a  valid  modulo  pQ   operation  requires  too 
many  details  which  it  is  undesirable  to  introduce  at  this 
point .   However ,  we  may  observe  that  the  circuitry  required 
is  minimized  if  m  and  consequently  pQ  has  the  least 
arithmetic  weight  possible  where  the  arithmetic  weight  is 
defined  as  the  number  of  non-zero  digits  in  the  canonical 
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(sometimes  called  the  non-adjacent )  binary  receding  of  the 

number o   Primes  of  the  form  2  -  m  with  m  as  small  as 

possible  and  of  least  arithmetic  weight  are  tabulated  in 
appendix  B  for  n  =  3  to  n  -  2*K 

Though  addition  and  multiplication  in  modular 
arithmetic  can  be  relatively  simple  and  fasts  the  analog 
of  division  in  integer  arithmetic  is  less  trivial  since  the 
process  requires  finding  the  inverse  in  FQ  of  the  divisor 
and  multiplying  it  with  the  dividend 0   Division  and  inver- 
sion do  not  have  the  same  importance  as  addition  and  multi- 
plication in  finite  field  operations 9  but  they  are  used 
frequently  enough  to  justify  studying  the  optimization  of 
the  relatively  complex  inversion  operation0   Since  a  search 
of  the  literature  indicated  that  only  table  look-up  has  been 
employed  to  any  extent  for  inversion,  it  seems  appropriate 
that  the  following  theorem  be  presented? 


Theorems   The  following  three  methods  are  valid  algorithms 

<c  p  01> 

for  the  computation  of  the  inverse  of  a  €  FQ, 


lo   (Product  Method) 

0 
Define  b^  $±      such  that  h±   €  {0,1},  3Q  -  a0  3±+1  =  &±   ©  (^ 

Pi         -i    n  0   b. 

and  pQ  -  2  -  )  ^2  0   Then  a  x  =  tt  ©  (^  )° 

i=0  X 
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2o   (Recursive  Method) 

If  a  =  1,  a"  =  1.   If  a  ^  1,   define  rQ  =  a  and  determine 

ns  q.   from  the  equation  p  =  q.  •  r.  -.  +  r.   where  n  is 

0 
such  that  rn  =  1.   Define  0  =  1,  Qi  =  (~qi+1)  ©  Qi+1  for 

0  <  i  <  n  -  lo   Then  a"1  =  QQO 

3«   (Factor  Synthesis  Method) 
If  a  ^  1,   calculate  q,,  r. ,  Q.   by  the  iteration 


p  =  qQ  -  a  +  rQ;   QQ  =  rQ  +  1 

0     .       0 
=  q±    •  ri_1  +  r±i      Qi  =  (Q1_1  $  (-1)  )  ©  q±9      i  >  0 


Find  n  such  that  r  •,  =  0o   Then 

a"1  =  (-l)n+1  0^9  (~l)n+1]  0  (rn)-1 

Proofs 

Part  1  is  based  on  Fermat's  theorem  which  states 

that  for  any  non-zero  element ,   a„   in  a  finite  field  of 

P0-l 
order  pQ»  a     =  L   This  implies  that  the  inverse  of  a 

P0~2 
is  a    o   The  proofs  for  the  other  two  methods  depend  on 

Induction  applied  in  a  direct  manner  and  so  are  not  included 

The  most  interesting  of  the  three  methods  of  in- 
version seems  to  be  the  recursive  method  since  with  it  there 
appears  to  be  the  ability  of  optimizing  the  procedure  by 
choosing  an  appropriate  pQ  and  it  also  seems  to  have  the 
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property  of  requiring  fewer  operations  than  the  other  two  in- 
version methods r  particularly  if  the  inversion  circuitry  is 
tuned  to  operations  in  PQ ,   Attempts  were  made  to  find  an 
analytic  result  for  the  maximum  and  average  recursion  depth 
for  a  particular  pQo   This  problem  appeared  to  be  untractable 
analytically?  so  the  brute  force  method  in  which  all  non-zero 
elements  in  the  prime  fields  generated  by  selected  pQ   of 
the  form  2  •=  m  where  n  ranged  from  3  to  16  were  inverted 
to  find  the  maximum  and  average  recursion  depths  for  these 
primes «,   The  results  of  these  calculations  again  did  not 
indicate  an  asymptotic  value  which  either  of  the  recursion 
depths  approaches o   There  does  seem  to  be  a  correspondence 
between  the  maximum  recursion  depth  and   aQ„   the  number  of 
factors  in  pQ  -  1.   The  calculations  seem  to  indicate  that 
the  more  factors  pQ  =  1  has,  the  smaller  the  maximum  re- 
cursion  depth  can  be  expected  to  be3   In  appendix  B  the  maxi- 
mum and  average  recursion  depths  are  listed  for  some  of  the 
tabulated  p  's. 

One  additional  property  of  the  prime  subfield 

which  might  lead  to  some  interesting  results  concerns  the 

2 
values  of  pQ  for  which  xt   *  1  is  irreducible..   If  this  is 

2 
the  case ,   x,  ■*■  1  will  be  a  prime  polynomial  and  the  vari- 
able x-,   in  the  finite  field  F,   will  have  the  property  of 
the  complex  variable  i0  This  points  toward  the  representa- 
tion of  models  involving  complex  numbers  in  which  i  is 
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2 
represented  by  x-,  o   If  the  polynomial  xr  +  1  is  irreduci' 

ble„   p0  -  1  is  referred  to  as  a  quadratic  non-residue 

(QNR)  and  if  it  is  reducible  to  two  linear  polynomials ,  a 

(12) 

quadratic  residue  (QR) o      This  is  the  reason  for  pro- 
viding two  primes  for  each  value  of  n  in  appendix  B<, 


If  a  significant  portion  of  a  computer's  time  is 
to  be  devoted  to  polynomial  manipulation  it  is  important 
that  the  operations  modulo  pQ  be  optimized .   Of  all  the 
operations  in  various  model  domains  these  prime  field 
operations  will  be  performed  most  frequently «,   Consequently 
it  is  expected  that  the  speed  of  the  polynomial  manipulator 
will  be  strongly  related  to  the  speed  of  these  prime  field 
operations  when  finite  field  based  modeling  is  used, 

5o3  plf  1  >  0 

We  have  already  seen  examples  of  p^'s  with  i  >  0o 
Specif ically,   (PqpP-.  =  x-,  -  l,p?  =  x2  -  l,„o.)   defines  an 
elementary  prime  context  sequence  in  which  each  of  the  fields 
F.   is  identically  equal  to  the  prime  subfield.   Although  the 
above  is  a  very  degenerate  form  of  prime  context  sequence,,  it 
is  expected  to  be  of  great  practical  value 9    the  reasons  being 
that  it  contains  a  totally  symmetric  set  of  polynomials.,   If 
the  modeling  is  performed  on  a  polynomial  in  several  vari- 
ables and  a  permutaion  of  the  variables  in  both  the  modeled 
and  the  model  polynomials  is  performed ,  the  model  resulting 


95 


from  the  permutation  will  be  identical  to  that  obtained  by- 
applying  the  modeling  function  to  the  permuted  modeled 
polynomial o   In  this  sense  the  modeling  is  not  affected  by 
the  permutation  of  the  variables  and  of  course  the  new  prime 
polynomial  context  sequence  will  be  identical  to  the  one  used 
prior  to  the  permutation0   Since  the  permutation  of  variables 
may  possibly  lead  to  procedures  for  optimizing  some  polyno- 
mial operations  this  symmetry  property  of  a  prime  context 
sequence  might  be  worth  persuing  further  for  other  less 
degenerate  sequences 0   It  is  also  Interesting  to  note  that 
the  above  elementary  prime  context  sequence  produces  re- 
dundant information  which  is  essentially  identical  to  that 

(8 ) 
produced  by  the  function  coding  used  by  Mart  in  <> 

Finding  prime  context  sequences  with  polynomials 
which  contain  more  than  a  few  terms  can  become  quite  com- 
plicated o  Fortunately  the  finite  field  theory  which  has 
been  developed  provides  us  with  one  main  theorem  which  can 
be  used  to  calculate  non= trivial  polynomials  for  the  sequences « 
The  theorem  is  stated  In  reference  (1)  theorem  19 „  pQ  135 » 
what  follows  is  a  repnrasal  of  that  theorem  using  the 
terminology  of  the  previous  chapters „ 


Definitions   Let  q±   a  pQ  °  X  2°    i  where  dQ  =  1  and 


<5o3o0i> 


d„   ■  degjpjo      Then     a   €  P-D    a  £  0     is   said  to 
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qi-l 
be  a  primitive  element  of  F.   if  a     =1  and  there  does 

not  exist  a  t   such  that  t  divides   (q,-l)   and  a  =  1« 


Theorem;   Let  a  be  a  primitive  element  of  F.   and  let  t 
-)o->°      *      be  an  integer  such  that  t  divides   (q.-l).   Let 
r  be  another  integer  such  that  g0c0do   (r„t)  =  lc   Define 

g  s  goCodo  (r9q,»l)o   If  4  divides  t  we  require  that 

t      r 
q.  =  1  (mod  h)  *     Under  these  conditions  i.  ,  -  a   is  an 

t  (Qi  =-1) 
irreducible  polynomial  belonging  to  the  exponent  — : — ± — -  0 

S 


The  properties  of  polynomials  used  in  the  prime 
context  sequence  for  practical  polynomial  manipulation  re- 
quire that  multiplication  in  the  finite  fields  be  optimized.. 
This  in  turn  dictates  that  the  prime  polynomials  contain  as 
few  non-zero  terms  as  possible »   The  binomials  produced  by 
the  above  theorem  are  therefore  optimum  in  this  respect., 

There  are  quite  a  number  of  open  questions  related 

to  finding  optimum  prime  context  sequences «   As  an  example  we 

9 
can  have  one  sequence  in  which    pQ  =  2  -9»P-.  =  x-,-  1  and 

3    2    3 
another  in  which  pQ  =  7,p,  =  x^  -  5  =  x^  +  3°   The  elements 

of  F-.   of  both  sequences  contain  approximately  9  bits  of 
informations  yet  the  structure  of  the  fields  is  quite  differ- 
ent 0   The  speed  of  operations  performed  in  F,   depends  very 
heavily  on  the  choice  made»   For  example „  if  operations  in 
Fq  require  negligible  time  compared  to  polynomial  operations 
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in  F..  0   the  first  sequence  should  be  chosen0   The  relative 
merits  depend  in  general  on  the  particular  methods  used  to 
implement  the  polynomial  manipulator . 
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CHAPTER  VI   CONCLUDING  REMARKS 

The  objectives  of  this  investigation  were  to  find  a 
class  of  redundancies  of  polynomials  in  several  variables 
over  the  integral  domain  of  integers  which  would  be  useful 
in  enhancing  the  power  of  polynomial  manipulators  by  in- 
creasing their  efficiency  or  their  capability  of  handling 
more  complex  problems  than  are  currently  possible „   The  finite 
field  based  modeling  functions  of  <3„6o07>  and  <*K2o08>  have 
been  shown  capable  of  doing  botho   The  most  illustrative 
example  presented  is  the  application  of  the  augmented 
modeling  method  to  the  determination  of  the  greatest  common 
divisor  of  two  polynomials?  it  is  shown  that  the  method  can 
be  used  to  increase  the  speed  of  determining  whether  two 
polynomials  are  relatively  prime  by  possibly  orders  of 
magnitude o   At  the  same  time  the  method  makes  feasible  the 
application  of  the  g0c0d0  test  to  polynomials  of  a  size  for 
which  it  would  be  impractical  to  apply  even  the  best  avail- 
able goCodo  method  without  using  modeling „   The  reason  for 
this  large  decrease  in  computational  effort  is  that  the 
coefficients  in  the  model  polynomials  are  finite  field  ele- 
ments which  implies  that  the  amount  of  storage  and  manipula- 
tive effort  required  to  work  with  them  has  a  relatively  low 
upper  bound o   Normally  the  problem  in  polynomial  manipula- 
tion is  that  the  size  of  the  coefficients  is  not  bounded  and 
consequently  the  storage  and  manipulative  effort  required  can 
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rapidly  reach  values  which  make  it  impractical  to  perform  the 
manipulation  even  on  a  fast  computer0 

Though  the  modeling  can  be  defined  so  as  to  produce 
conclusive  results  in  most  practical  cases „  it  is  to  be  under- 
stood that  there  will  always  be  cases  in  which  the  results  are 
not  conclusive o  a  fact  which  must  be  borne  in  mind  when 
polynomial  manipulations  dependent  on  decisions  made  from  the 
models  are  used.  Although  much  of  the  information  in  the 
polynomial  can  be  retained  in  its  model 0  the  possibility  of 
decisions  being  based  on  lost  information  does  exist  and 
hence  the  models  cannot  be  manipulated  blindly  without  taking 
this  into  account „   The  virtual  degree  augmentation  to  the 
modeling  is  one  example  of  a  technique  which  may  be  used  to 
prevent  this  type  of  decision  errorQ   Inconclusive  results 
may  still  be  encountered,,  but  with  the  augmentation  they  are 
detectable,  making  it  possible  to  take  some  remedial  actions 
such  as  reordering  the  variables „  using  a  different  modeling 
or  even  changing  the  prime  context  sequence o   If  none  of 
these  actions  produce  a  conclusive  result  one  can  always 
return  to  the  modeled  domain  and  apply  the  algorithm  without 
using  modeling  provided  of  course  that  this  is  a  practical 
possibility o 

It  is  postulated  that  the  same  techniques  as  used 
in  the  g0c0d0  calculation  can  also  be  successfully  utilized 
to  reduce  the  computational  complexity  in  elimination  and  in 
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the  reduction  of  systems  of  differential  equations;  results 
in  this  area  are  expected  to  be  particularly  worthwhile  in 
regard  to  making  algebraic  manipulation  efficient  enough  to 
be  of  more  general  use  than  is  now  the  case0   Some  additional 
areas  in  which  modeling  might  prove  useful  are  polynomial 
factoring9  error  checking,,  operation  optimization  by  model 
determined  orderings  and  printed  polynomial  output  monitoring c 
There  is  also  the  possibility  that  the  modeling  can  be  applied 

to  functions  other  than  polynomials  in  a  manner  analogous  to 

(8 ) 
the  coding  used  by  MartinD 
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On  the  following  page  a  polynomial  in  I~  with 
64  terms  of  one  decimal  digit  each  is  mapped  Into  all  its 

possible  model  domains  under  the  prime  context  sequence 

2 

P0  =  70  P-,  =  x^lfi  p2  =  Xp-10  p~  =  x~-lo   The  information  is 

compressed  in  such  a  way  that  the  polynomial 


x-l  +    (2+3x^)x2  +   ( 4x^+5x^2+6x2  )x3  +    (7x1+(8+9x^)x2)x^ 


would  be  represented  by 
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An  X  in  a  position  indicates  that  the  coefficient  which 
should  be  present  is  greater  than  9o 
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*   indicates  the  start  of  a  polynomial  in 
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APPENDIX  B 


Minimum  Weight  Primes 
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=  arithmetic  weight  of  pQ 


oQ  ss  number  of  factors  in  P0~l 

MR  =  maximum  recursion  depth  in  recursive  inversion  method 

AVR  -   average  recursion  depth 
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